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ABSTRACT 

Parameter plane techniques were first introduced in an IEEE paper 
dated November 4, 1964. The paper dealt mainly with the theory of the 
parameter plane whereby the roots of a pottynoditel could be determined 
graphically in terms of two parameters which may appear linearly in 
any of the coefficients. 

In this text, parameter plane techniques are applied to the 
compensation of linear feedback control systems by both graphical and 
analytical means. Parameter plane equations are extended to include 
parameter products and three parameters. An attempt is made to show 
the complementary roles of the parameter plane and root locus. 

The writer wishes to express his appreciation for the assistance 
and guidance given him by Dr. G. J. Thaler of the U. S. Naval Post- 


graduate School in this investigation. 
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ds Introduction. 

The analysis and synthesis of feedback control systems, or the 
compensation of same, can be effected by three general methods. The 
first of these can be called the integral criteria. Here a cost function, 
in which is inherent the system design specifications, is minimized with 
respect to certain variable system parameters. This method is mainly 
applied to the statistical properties of feedback control systems. The 
second method is the Bode frequency response method whereby the systems 
open loop transfer function is manipulated to obtain the desired system 
response. This method has its inherent weaknesses, such as difficulty of 
application to non-unity feedback control systems, difficulty in inter- 
preting closed loop transient response in terms of open loop frequency 
response, difficulty of varying more than one parameter, and the short- 
comings inherent in the approximations which are required when applying 
the frequency response method to compensation problems. Third are the 
algebraic methods. Under this heading can be listed the root locus method. 
The shortcomings of this widely known and widely used method are familiar 
to its users. Its greatest disadvantages lie in plotting the actual locus 
of Gracy and in the fact that only one parameter can be varied conven- 
iently. In references (1) and (2), Ross-Warren and Pollak respectively 
developed algebraic methods of cascade compensation using root locus 
techniques, but the inherent disadvantages of the root locus were still 


present. In reference (3) Mitrovic developed an algebraic and graphical 


aa computer program is introduced in section (5) which makes this 
disadvantage less significant. 








method of obtaining the roots of a polynomial in terms of two variable 
parameters. Later in reference (4) Ohta developed some additional 
sketching techniques which greatly facilitated the plotting of the 
Mitrovic curves. In references (5) and (6) Choe and Hyon respectively 
applied and extended the Mitrovic method to the compensation of linear 
continuous feedback control systems. The inherent disadvantage of the 
Mitrovic method is that the variable parameters may appear in no more than 
two coefficients of the characteristic equation, which reduces the flex- 
ibility of the method. In reference (7) Siljak introduced a method of 
obtaining the roots of a polynomial in terms of two variable parameters 
which can appear in any and all the coefficients of the polynomial. In 
this text, Siljak's method is applied and extended to the compensation 
of linear continuous feedback control systems. General methods of com- 
pensation will be developed and an attempt will be made to relate the 
root locus, and the parameter plane as a set of complementary techniques 
which when applied in conjunction with one another represent the most 
adequate tool to date for solving the problem of compensation of linear 
feedback control systems. 

The relationship between being able to place the roots of a polynom- 
ial at specified locations in the S-plane and the compensation of feed- 
back control systems is as follows. The basic idea is that any feedback 
control system, including any added compensators which may contain vari- 
ables, can be reduced to or can be represented by a ratio of two poly- 
nomials which is the closed loop transfer function. Well known methods 
are available whereby a specified system response, in terms of over- 


shoot, bandwidth, settling time, steady state accuracy, etc., can be 








obtained by placing a pair of complex conjugate roots of the character- 


istic equation at a specified location in the S-plane, while ensuring that 
the real part of this complex root pair (called the dominant roots) is 
smaller in magnitude than the real parts of the remaining roots of the 
characteristic equation. The problem of compensation therefore reduces 
itself to one of moving the roots of the characteristic equation to the 
desirable locations. The usefulness of the parameter plane to achieve 


this will soon become apparent. 








oe Derivation or the basic parameter plane equations. 
A feedback control system's characteristic equation can be re- 


presented as a polynomial of the following form: 


yy 


f(s)= £¢ aS = 0 (ai) 
K=O 
Where the coefficients a, (k = 0, 1,...-,m) are real, and S is the 
complex variable S = -G+ jw = - y w+ jw 1- 5 ; (2-2) 


w is the undamped natural frequency, and 7wis the relative damping 
coefficient. It is noted in reference (7) that sk may be represented 


by the following: 


ss awit (- 7) +4 f1- 7* we 7) (2-3) 
k k 


where T.(- 7) = CIC F) and (FI = CD UF). ACF) 
and UC 7) are Chebishev functions of the first and second kind respect- 
ively. Values of zeta will be considered such that 0 = i = 1 and 
values of w such that 0 = w <s o . The values of Tk and UL are tabu- 
lated in Appendix I for various values of zeta. But more useful to digi- 
tal computer employment, they can be obtained from the following recursion 


relations: 


se lem 7 ) * Ty. 1a 2) 


U4 07) -2U07) +0, 07) = 0 


Here TC 7) -1,T,07) = 7, Uo 7) = 9, UC FI =1 
Substituting equation (2-3) into (2-1) and setting the real and imagin- 


ary parts to zero independently one obtains: 
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vy} 


k 
2a Dales) 7 nse 


y4 k 
Z aU (- J) = 0 


K =< 


(2-6) 


Employing equations (2-5) one obtains from equations (2-6); 


Wr 


2 (-1) ‘awit, _5( py Sw 


fi 
= 
K =. 


=e) 


(-1)"a wt, (  ) = 0 


(2-7) 


Consider the coefficients ay. of the characteristic equation (2-1) as 


linear functions of the variable system parameters as follows: 


Seep ee 


(2-8) 


Employing the above relation for ays equations (2-7) give the following 


relations: 
/ aa i— = 
Where: 
- 2 k, ok ei AE 
“17 a Soe ees Bee Okman 
m M 
E k _ 4k, k ‘ 
C= z,¢ 1)"e, UL, C, eX 1) “c.w UL (2-10) 
k, k ep k, k 
= = = > -_ 
D, ZS 1) dw kel Do =, 1) dw UL. 


Since equations (2-9) are two linear equations in 


and beta, Cramer's rule can be applied to obtain: 


C,D,-CoD, 
ByC,-B,C, 


ol = Sgr lee 


122. 1 
Equations (2-11) are now functions of zeta and w. 


varying zeta or by fixing zeta and varying w, the 


yi 


B,C,-B,C 


the two unknowns alpha 


(2-11) 


Hence by fixing w and 


constant w or constant 








zeta 5S plane contours respectively can be mapped into the real domain 
of the alpha beta plane or parameter plane. 


In reference (7) the following relationships are utilized: 


k 2 
Ss = Fae he a. 
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PL hp + 2w Pe + ow Peel = 0 (2-12) 
9 + 2wQ, + ae = 0 
et 1 i woke 
P=, Ey eZ meer 0; Qt 
P =-w Q -w°0 
k / k k-1 
Here Pe and % are related to the Chebishev functions by 
sage _ k k 
PrteMat CooFuetC-1) lage) (2-13) 


Rind eke Fidun Cy Rthyerhy 3) 


By using equations (2-12), and (2-13), one obtains proceeding as before: 
m m 
ane Opey He €.*en ° e-®) 
Employing equations (2-8), (2-14), along with Cramer's rule one again 


obtains equations (2-11) where the following expressions now apply: 


an m 

~ 2 ee B= Zee 
Vi m 

cited Se Cale 2a, =, %K% a 
m mM 

ea Do ES 





Beuations (2-11) and (2-15) are useful for mapping constant zeta-omega 
curves from the S-plane into the parameter plane. As will be seen later 
these curves play an important role in dominance considerations. 
If the complex variable S is substituted in equation (2-1) by lett- 
ing 
S=-6 (2-16) 
where sigma corresponds to values of S§ along the real axis, then in 
accordance with equations (2-8), the characteristic equation (2-1) be- 
comes: 
Le (oy e*+ BF cst + = (-1) “a, = 0 (2-17) 
K=O K=0 K=0 
The above represents a straight line in the alpha-beta plane for a given 
value of sigma. Hence a point on the real axis in the S-plane maps into 
a straight line in the alpha-beta plane. Also for a given value of alpha, 
beta, and sigma which satisfies equation (2-17), then the characteristic 
equation (2-1) must have a real root at minus sigma. For the constant 
zeta and omega curves as defined previously, if for certain values of 
alpha and beta, say for a value obtained from equations (2-11) with a 
certain value of zeta and omega, then the characteristic equation (2-1) 
has a pair of complex roots at S = - 7 w + iv] 1- ji 
It is important to note that by applying equations (2-11) and (2-17) 
one can, for a specified value of zeta, omega, and sigma, compute the 


value of alpha and beta so that the characteristic equation will have a 


as 
pair of complex roots at say S = ~ 711 = Jw, |. 71 » and a real 
root at S = - 61° The m-3 remaining roots of the characteristic equa- 


tion can then be determined by dividing out the three known or specified 
roots. This method where zeta, omega, and sigma, or just zeta and omega 
are specified, and the computations for alpha and beta are done algebrai- 


cally, will be referred to as the albegraic parameter plane solution. 


6 








To solve the problem in general, for all values of zeta, omega, and sigma, 
it is necessary to plot a family of parameter plane curves for various 
values of zeta, omega, sigma, and if desired, zeta-omega. On the result- 
ing parameter plane plot one can, by picking an M point or operating 
point, graphically read from the curves the values of the m roots corres- 
ponding to the mee order characteristic equation. This latter method will 
be called the graphical parameter plane solution. 

The algebraic solution has the advantage that the labor of plotting 
the curves can be evaiaedme but it has the disadvantage that without the 
curves it is sometimes difficult to pick the most optimum value of zeta 
and omega so as to ensure dominance and still meet the system specifica- 
tions. The graphical solution has the advantage that one has a picture 
of the way the roots of the characteristic equation move around in the 
S-plane as alpha and beta are varied. This enables one to pick the values 
of alpha and beta corresponding to the best values of zeta, omega, sigma, 
and zeta-omega for all the roots of the characteristic equation. This 
latter feature of the parameter plane points out a strong argument for 
trying to obtain the parameter plane curves. If a digital computer is 
not available then by using the relationships derived in section (3-3) 
under sketching techniques, along with a desk calculator or slide rule, 
the curves can be plotted with some labor. Under these circumstances it 


is questionable whether the algebraic or the graphical solution would be 


better. Which one is used is a matter of personal preference. 


othe computer program presented in section (6-1) is most helpful in 
reducing this labor. © 
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Dig Application of the parameter plane to the compensation of linear 
continuous systems, 
3-1. Algebraic solution, 

In this section it will be assumed that the system performance specifi- 
cations have been given in terms of placing a pair of complex roots at a 
specific value of zeta and omega, say OO and Wy> with the error coef- 
ficient K. being greater than or equal to a specified value. If the specifi- 
ed location of the roots is such that after computation of the necessary 
value of alpha and beta, the remaining roots of the characteristic equa- 
tion are located so that the specified roots are not a dominant pair, then 
either a different value of zeta and omega will have to be used (possibly 
at the sacrifice of some measure of the system performance), or a differ- 
ent method of compensation will have to be used. 

In section (3-2) a method is presented whereby the dominancy specifi- 
cation can be met by introducing a third parameter. 

3-1-1. Feedback compensation. 

Figure (3-1) represents a unity feedback control system. In order 

to meet the system specifications a feedback compensator H will have to 


be used. Let 
G=kK = K (3-1) 


Where K is the forward path gain which can be varied and e(S) is a poly- 
nomial in S representing the poles of the open loop transfer function of 
the uncompensated system, The letter L in equation (3-1) corresponds to 
the system type. For a type O system, L = 0, for type 1, L= 1, for type 


2, L = 2, etc. By definition, the error coefficient is given as follows: 
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. re we Th Oe oma 






a = 


— einen 











7. Stee 





, 


| 





eee. e (3-2) 
e ce 


S70 
Here Gon is the open loop transfer function of the compensated system. 
Sometimes, as for example in reference (8) the error coefficient is design- 
ated as — for a type zero system, K, for a type one system, and K. for a 
type two system. 


Tachometer plus acceleration feedback. 


Here one lets 
Z 
H = KS oF KS (3-3) 
The resulting compensated system's characteristic equation becomes: 
e(S)4+K+(K,SHK 8°) = 0 (3-4) 


or by expanding e(S), equation (3-4) becomes: 


m m-l 2 
S + e 1-15 +...+(e,+KK,)S +(e, +KK, )Ste +k = 0 (3-5) 


where L is taken to be zero for a type zero system which one can consider 
as the most general case. The following results also apply to a type 
one system if ey is set to zero, or to a type two system if both ey and 


e, are set to zero, etc. From equations (3-2) and (3-4) the error co- 


1 


efficient becomes: 





K = ain s° K = K (3-6) 
2 
S40 e(S)+K(K, S4+K,S2) ey 
or if the uncompensated system is type one: 
Ke = K @=7) 
e, TKK. 
and if the uncompensated system is type two: 
re K (38) 
ie 


Note: If the uncompensated system is type two, the compensated system 


10 





would be type one if tachometer feedback or tachometer plus acceleration 
feedback is used. 

In the compensated system's characteristic equation (3-5) let alpha 
= KK, and beta = KK Equation (3-5) then becomes: 

sy er et ieee Cee +B Ste +K=0 (3-9) 

m-l1 2 1 oO 

Recalling from equation (2-8) that in general the coefficients of the 
characteristic equation are of the for: 


k 


equation (3-9), m = k, one finds that: 


Al = bl “ c¢. B +d, and corresponding] y. th 


d = e5 + k, bo =c =0, d, =e,, b, = 0, c, = 1, d, =e, b, = 1, 


oO O i 1 et 1 2 2” 92 
oo OF aa a dw Dey = 0, C4 = O, etc. 
Then from equations (2-10) one obtains: 
. Due 5 A ee 
By = (-1) w U, =w B, W U, 
Cy = “WU 5 = 0 Cy = “wu, = -w (3-10) 
ms k, k kk 

ay =, Earn ea a z, (-1) aw UL, 

Use was made of the fact that U_, = -l, U, = 0, and UL, = 1 (see appen- 


dix I-B). Using the expressions for alpha and beta as given in cquations 


(2-11) along with the above information the following relations evolve: 


Y! 
: ee YS 
= c,d, C.D, _ WwW ek 1) dw 


B,Co-BoC, We 


k m 
k=) =) eee (ay we 2y (3-11) 
K=0 k k-1 


m 

| B = =. (-1) “aw *(U,-U,U,_4) 

At this point alpha and beta may be linear functions of K, the for- 
ward path gain, and one can use the steady state error specification to 
put K in terms of alpha and or beta. Since zeta and omega were assumed 
to be specified, then equations (3-11) can be used to solve for alpha and 


en 
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beta. Then K and Kh. can Treadisyebesdcter mained, 


Example 3-1 
The systen given in figure (3-2) is to be compensated using tacho- 
meter plus acceleration feedback. The system specifications are as 
follows: 
1. Complex roots at 7 = .7, w= 10. 
2. Ky 2 6, not to be reduced. 


Solution: 


From equation (3-2): 


K_ = lim sig = K >6 
S>o0 a 2+ RK 


Oe K —12 + 6KK, . The compensated system's characteristic equation is: 


s> + $23 4 KK.) + S(2 + KK.) + K = 0 (3-12) 


Letting alpha = KK, and beta = KK equation (3-12) becomes: 


MeEeaGERe ICE x = 0 (3-12a) 


Employing equations (3-10) the following can be obtained: 


By = 100 Cc) = 0 
dD, = -1100-K B. = 140 
C, = -10 dD, = -1120 


Using equations(2-11) one obtains: 
o& = 10(-1100-K) 8 = 140(-1100-K) + 56000 (3-13) 
- 1000 -1000 
Note: The preceding expressions could have been arrived at directly by 
employing equations (3-11). From the steady state accuracy specifica- 
tions it is necessary that: 


m. =12 +6 6 , hence let K = 12468 (3-14) 
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Using equation (3-14) and (3-13) beta is found to be: 


C 


= 
: 
ot 


140(-1100-12-6 © ) + 56000 = 625 
-1000 


Therefore K = 12 + 6(625) = 3762 and oC = .01(1100 + 3762) = 48.6 


Since cm = KK. then K = 48,6 = .0129 
“3762 


Also from 6 = KK, it is seen that = ,166. 


The compensated system's characteristic equation becomes: 


s? 4 51.628° + 627S + 3762 = 0 (3-15) 


Now zeta = ./ and omega = 10 corresponds to 5? + 14S + 100. When this 
quadratic is divided out of equation (3-14), the remainder is S + 37.62. 
Hence zeta-omega of the desired roots = 7 <<37.62, and the complex roots 
are dominant so the problem is solved. 
Example 3-2 

The problem presented in example (3-1) will now be solved by intro- 
ducing the error specifications at the beginning of the solution instead 
of at the end. Equation (3-12) is as follows: 


s° 4 §°(3 + RK) + S(2 + KK) + K = 0 


Again let alpha = KK, and beta = KK Then from the steady state error 
requirement one obtains: 
K=12+6f (3-16) 
Substituting this expression for K into equation 3-12 results in: 
so +5773 +00) +5(2+8) +124 68 =0 


Therefore bo = 0, c= 6, do = 12, and all other coefficients are as 


before. Hence from equations (3-10) it is seen that: 


B, = 100, C, = -6, D, = -1112, B, = 140, C, = -10, D, = -560 
and of = -7760 = 48.6 p = 625 K= 3762 
-160 
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The above result agrees with example (3-1). 
Note; Equations (3-11) could not be used here since they are based on 
applying the accuracy specifications at the end. 
Tachometer feedback only. 
Let H = KS (3-17) 
The characteristic equation of the compensated system becomes: 
Se $e,8 +(e, +) Ste ie (3-20) 


Oo 
From equations (2-10) one obtains: 


By = Be = =—W 
c= -1 Gc, = 0 (3-21) 
p= & (-1)*awt D, = = (ap aan 

1 K=O k k-l 2 Ke= k k 


From equations (2-11) it is seen that: 


ma 


oh = = (-1)“aw'u, and SB = Z, (-1)Saw't, , (3-22) 

If a given zeta and omega are specified, the alpha and beta can be 
computed from equations (3-22). The error coefficient is then determined 
directly from equations (3-6), (3-7), or (3-8). Hence the error coeffi- 
cient is fixed once zeta and omega have been chosen, so if a certain 
error specification is to be met, the specified values of zeta and omega 
may have to be adjusted to meet it. 

If the error specification was the overriding specification to be 
met, then zeta could be fixed at some reasonable value. By means of the 
given Ka» alpha could be computed from equations (3-%), (3-7), or (3-8). 
Equations (3-22) could then be used to solve for first omega and then 
beta. The calculations would be more tedious however. 

Example 3-3. 


Figure (3-3) shows the same system as used in the previous two ex- 


amples only now tachometer feedback alone will be tried. The same system 
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specifications are to be met, i.e., K, Zz 6, zeta = .7, and omega = 10. 
The compensated system's characteristic equation then becomes: 


si 4 357 + (2 + KK,)S + K = 0 (3-23) 


Using equations (3-19) and (3-23) one obtains: 


s° +357 + (24% )S +P =0 (3-24) 


From equations (3-22) alpha is found to be: 


Oh = -2 + 30(1.4)-100(.96) = -56 

Since alpha is negative it is seen that positive tachometer feedback is 
required. Since the coefficient of the first power of S in the character- 
istic equation would then be negative, the system would be unstable. 
Hence the desired system specifications can not be met with tachometer 
feedback alone. 
Acceleration feedback only. 

Let H = KS" (3-25) 
The characteristic equation of the compensated system then becomes: 


a eb St (e, KK, )s* +e 


ae Ste +K=0 (3-26) 


1 
Let alpha = KK, and beta = K. (3-27) 


Then from equations (3-26) one obtains: 


m m-1 2 7 
S + ewe? nee (e, +o )S” + e,S + ey +6 = 0 (3-28) 
Using equations (2-10) and (3-28) results in: 
ae Ee A ee 

By = Ww Us Ww B. = Ww U, 

C, = U1 ed C. = 0 (3-29) 
yy k, k rn k 

D, = =f» dw Uy dD, ee dw'U, 


Solving for alpha and beta results in: 
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m k, k-2 
sills i 7 ] S (-1) dw UL. (3-30) 
- U Kea 
WU, 2 
k. ok m ee ks 
GP = 2,60 dwt, 7h = (-l) dw, 
= Ue K=O 


Calculations for alpha, beta, and the error coefficient are performed 
in the same manner as with the preceding tachometer feedback calcula- 
tions. 
Example 3-4 
The system of examples (1) and (2) will now be compensated using 
acceleration feedback as indicated in figure (3-4). As before K, = 56. 
zeta = ./, omega = 10. Therefore K. = K/2 and the error coefficient is 
unaffected by the acceleration feedback. Hence one can choose K = 12 to 
meet the specifications. The compensated system's characteristic equa- 
tion then becomes: 
Sa + (3 + KK,)S* +25 +K =0 (3-31) 
If K in equation (3-31) is set to its prescribed value of 12, only 
one parameter remains and the parameter plane equations produce an indeter- 
minate solution. Therefore K will be left as the variable beta. Since 
beta is fixed by the chosen values of zeta and omega, then so is the error 
coeificient, and it will most likely not agree with the error specifica- 


tion. In view of this, the solution proceeds as follows. Making the 


usual change of variables in equation (3-31) results in: 
3 2 
S~ + (3 +00 )S~ + 25 a =10 (3-32) 


By employing equations (3-30) one can solve for alpha and beta. 


of 


P 


-1_ [.-2/10 + 3(1.4) -10(.96)] = 4 


1.4 
3(10)2 -(10)7°(1.4) -1A4 [-2(10) + 3(10)7(1.4)-(10)7(.96)] 


or 8 = -700 
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Fron the negative value of beta it is concluded that the desired 
values of zeta and omega cannot be obtained using acceleration feedback, 
and of course neither can the desired error specification be obtained. 

One would therefore choose another method of compensation. 

If one chooses to use feedback compensation then perhaps tachometer 
plus acceleration feedback should be tried first using equations (3-11) 
and the appropriate steady state error specifications. If the specifica- 
tions cannot be met in this manner then it is obvious that they cannot be 
met by either tachometer or acceleration feedback separately. In this 
case either the system's specifications must be modified or another 
scneme of compensation must be employed. If it is found that the specifi- 
cations cen be met by combined tachometer and acceleration feedback, then 
if desired, equations (3-22) or (3-30) can be used to see if tachometer or 
acceleration feedback alone can be used. 

Cases where feedback is not around the forward path amplifier. 

Figure (3-5) illustrates a compensation situation that sometimes occurs 
in practice. This is the situation where it is not possible or practical 
to get at the input terminals of the error detector and the feedback has 
to be inserted at the output terminals of the amplifier. The solution to 
this problem is solved by means of an example. 

Example 3-5 

Figure (3-5) shows the system that was used in example (3-1) only 
now the feedback is inserted at the output terminals of the amplifier 
represented by the gain K. The same ees specifications are to be met 
ime. 5 K, = 6, zeta = .7, and omega = 10. The characteristic equation now 
becomes; 


ge aa K,)s" + (2 +K,)S +K = 0 (3-33) 
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Letting alpha = KA and beta = K, equation (3-33) becomes: 


s+ (3 +6 57 + (2+8)S+K = 0 (3-35) 
Comparing equation (3-35) to equation (3-12a) it is seen that they are 
identical, so the solution obtained for alpha and beta in example (3-1) 
applies. This points out an important advantage of the parameter plane 
method. This is that the solutions depend only on the characteristic 
equation and not on the . svstem that the characteristic equation was 
formed from. 

From example (3-1) it was found that alpha = 48.6, beta = 625, and K 
= 3762. In this example there are no additional computations necessary 
to find K. and Ke» since these are now the system's parameters alpha and 
beta. So K = 625 and K, = 48,6, 

This same general principal can be applied to control problems in- 
volving tachometer feedback alone or acceleration feedback alone. 

3-1-2 Cascade compensation. 

Fieure (3-6) represents a unity feedback control system which in 
order to meet the system's specifications a cascade compensator G. is re- 
quired. Let G have the form of equation (3-1). 

i.e.g G=K = K (3-1) 

ey s" +e olay, ce gh 
m-1 L 
K is the forward path gain which can be varied and e(S) is a polynomial 
is S representing the poles of the open loop transfer function of the un- 
compensated system. The letter L again corresponds to the system type. 
Let: 


G = P(S + Z) (3-36) 
P 


This compensator has a D.C. gain of unity so its placement in the forward 


path will not affect the steady state accuracy. It is assumed that the 
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uncompensated system's forward path gain had previously been adjusted to 
give the correct steady state accuracy. Using the form of G. as indicated, 
the values of Z and P are computed to give the desired system response. 
If P is less than Z% a lag network is needed and the factor P/Z of the 
compensator is inherently present due to the physical nature of the compen- 
sator, which is assumed to be an R-C network. See appendix II concerning 
the nature of lag-lead R-C networks. In this case all forward path ampli- 
fier gains can remain unchanged to meet the stipulated accuracy specifica- 
tions. If however, the computed values of Z and P are such that P is 
greater than Z, a lead network is required and the compensated system's 
forward path gain will have to be raised by a factoz of P/Z. The physical 
nature of the lead network is such that the factor P/Z is not inherently 
present, so to maintain steady state accuracy this factor will have to 
be provided either by adding an amplifier in cascade with the lead net- 
work or by raising the gain K of the existing amplifier by this factor. 
Continuing then, the procedure is as follows. The compensated system's 


forward path transfer function is: 


a™ GiG) = Raa PS +12) Se sirme en es (3-37) 
: a(sy Zz S +P S +P e(S) 
Applying the definition of the error coefficient to the compensated system 


one obtains: 


K = Lim f KK. ¥(S + P/¥ | = K_ (3-38) 
e 
L 





cm o8¢ e(S) (S + P) 


Again assuming a type zero system where L = 0, the compensated system's 


characteristic equation becomes: 


e(S) (S +P) +K ¥ (S + P/¥) = 0 (3-39) 


or in general after expanding equation (3-39): 
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m-l 


otrr l 2 o~ 
+ e 7-2» ieee (Pe, + e,)5 + (3-40) 


n 
+ (P + eae)? + (Pei 


(KY + Pe, + e)s + P(e. + K) =0 


1 
Letting alpha = P and beta = ¥ , equation (3-40) becomes: 


m-l1 


m+1 | m 2 
S + (A+ e-)8 (e 1% + e-9)8 +...4+(e, 4 + e,)s (3-41) 


+ + + = 

(KB e, A es A(e. + K) 0 

Comparing equation (3-41) to the general form of the characteristic equa- 
tion as specified in equations (2-1) and (2-8), it is apparent that K = 
m+ 1 (the order of the equation), bo oot K, c #=d 2=0, b, =e 


c, =K, d, =e , b, =e,, c = 0, d, =e,, etc. 


2* “9 2 1 
It is important to note that the parameter plane variable beta re- . 
presents the pole to zero ratio of the cascade compensator. In references 
(1) and (2), Ross-Warren and Pollak respectively, utilized the concept of 
root relocation zones to divide the S-plane into regions where lag compen- 
sation or lead compensation is required. By assigning variables in the 
above manner, the parameter plane is effectively divided into corresponding 
regions above and below the straight line beta = 1. For values of beta 
less than one, a lag network is required and for values of beta greater than 
one, a lead network is required. In addition if beta is greater than 10 
or less than .1, a multiple1ecdor multiple 1az ‘respectively is required. 
A multiple section enn will also be reived if the computed value 
of either P or Z turns out to be negative. A method is given in section 
(4-6) for the design of double section compensators. If more than two 
sections are required, the compensation has to be done in steps. 
In this case complex roots are placed at some intermediate value of 
zeta and omega and a new characteristic equation is then computed. An- 
other compensator can then be designed on the basis of this new character- 


istic equation to place the roots at the desired value. Thus by compensa- 
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tion in steps in conjunction with the double section design of section 
(4-6), a four section compensator could theoretically be designed. The 
use of more than four sections is questionable and in this event it would 
be better to employ combined cascade and feedback compensation, or per- 
haps feedback compensation alone. 


On the basis of equations (3-41) and (2-10) it is found that: 


Z , 2 _4\k-2 k- k-1 k-l 
Bp = (eg +R) twee, t..tIy Ww + CWE 
Cy = 0 
we K-20 k-2 kK-l kel, 4 
D, =we, +...+(-1) en" U3 + (-1) e-1" UL-9 
(-1) wt, (3-42) 
_ 2 k- k-1. k-1 
By = -we, +wW eoU, +..2+(-L) vies + (-1)° ow UL 
Cy = -wK 
he 2 k-2 k- 
dD, = “we | +w e,U, +e. «has eo" aU 
k-l] k-] k 
+(-1) el Und + (-1l) wD, 
and from equations (2-11) it is found that: 
cm? B= lee (3-43) 
ey eee 


In equations (3-42), one sets e, = 0 for a type one system, @ ar Cae 0 
for a type two system, etc. On the basis of equations (3-42) and (3-43) 
a cascade compensator can be designed. 
Example 3-6 
Problem: 

Design a cascade compensator for the system shown in figure (3-7) 
to place a pair of characteristic roots at zeta = .5 and omega = 1. The 


error coefficient K should be 50. 
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Solution: 

From figure (3-7) it is apparent that K, = K/10 = 50 or K = 500. The 
characteristic equation is: 
s* + S703 +P) + s*(17 + 8P) + S(10 + 172 + KY) + P(10 + K) =O (3-44) 
Letting alpha = P and beta = {,equation (3-44) becomes: 
s* + s°(8 +h) + $7(17 + 80L) 48(10 + 17% + 5008) + (10 + 500)% = 0(3-45) 


Applying equations (3-42) and (3-43) one obtains: 


BL = -503 Bo = -9 

Cc, = 0 Co = -500 

D, = 9 D, = 6 

f& = .0179 = P 6 = .0117 = ¥ 


But Y = P/Z so Z = 1.529. 


G. then becomes: -O117 (S -+ 1.529) 
(S + .0179) 


This is a lag network Divera the factor .0117 is inherent in the R-C filter 
design. Due to the small value for gamma, the size of the filter compon- 
ents may be unreasonable, and a double lag network could be designed us- 
ing the method of section (4-6-2). Instead the problem will be solved in 
section (3-1-3) using combination cascade plus tachometer feedback compen- 
sation which will permit the use of a single section lag network. 

3-1-3 Combination cascade and feedback compensation. 


Derivative signal not enclosing a cascade compensator. 


In figure (3-8), G = K/e(S) and G. = 8S + 
S + 


P (3-46) 
P 


First let H(S) = KS: Then in view of equations (3-46) the compensated 


system's forward path transfer function becomes: 


G ¥S+P . K/e(S) = YS+P. K (3-47) 
oT 1 + KH(S)/e(S) S$ +P e(S) + KH(S) 


Zi 
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and Ko = lim s& G = K (3-48) 


320 si e, (S) + Lim KH(S)/S™ 
S70 
Note: The system type will change if the lowest order of the deriva- 
tive signal fed back is less than the system type number. A few error co- 


efficients are given in table (3-1) for different types of uncompensated 


systems, In the table let L be the type number of the uncompensated 


system. —, 
Table of K 
e 
1 a 0 1 2 
H(S)J, 
KS K/e. K/(e,+KK, ) 1/K, 
2 
K,S K/e. K/e, K/(e,4kK,) 
Table (3-1) 
The compensated system's characteristic equation is: 
Se(S) + Pe(S) +SKH(S) + PKH(S) + K ¥ S + PK = 0 (3-49) 


Using equation (3-1) for e(S) and letting alpha = P, beta = ¥ , and 


k = mtl1, equation (3-49) becomes: 


sk - (XK + ee | + (e__y%+ om on +...+(e, + KK, + e, X)S* 
+ (e,K + KL + KG + e.)S tol (K +e.) = 0 (3-50) 


The parameter plane variables are then: 


Br ise-e -K + we, feet) ee wy ot (1) Py © tes 
Cc, = 0 
D, wate, + RK Dw ew Ut. t(-1) 76 wu 
“ (-1)"te wu + (-1) wd, (3-51) 


30 






a i Ana 





— 9 en aii ; 


= —_—= - ee ae ==> 


oa ina Cet ew Us te tD ew 7U 
" cay ty 

Cy = -Kw 

D, = -ew “- (e, + KK, )w U, tect) Fe wh ?u, 
+ (-1)"e wi tu is (-1)Sw"u, 


In terms of equations (3-51) one can solve for alpha and beta: 


oh = -D./B, B = Vizag (3-52) 


~WKB, 


One can now compare the above expressions for alpha and beta with 
equations (3-42) and (3-43) to see the effect of tachometer feedback. 
Let Ba» Ci» Dy» Bos Cos and Ds represent the quantities given by equations 
(3-42) where only cascade compensation was used. Then in terms of the 


primed quantities, equations (3-51) become: 


By By B, By KK LW 
C, C, 0 C, Co (3-53) 
D, = D! + KK a D, = Di + KK yall 

1 1 t Ze 2 t Z 


Let the expressions for alpha and beta as given by equations (3-43) 
be designated of ' and B ', Then in terms of these quantities, equa- 


tions (3-52) can be expressed as: 


a 2 ' = i 2 ’ - 
of = ( D) - KKYw )/B, of KK .w /B5 (3-54) 
Z 2 -R ? g 2 2 ’ g 
S = (-Bowk + KD) + Kikw )/B) + Kw, +2 


Now it can be seen from equations (3-54) how tachometer feedback modifies 
the values of “' and er as computed for cascade compensation alone. 


For instance if the pole to zero ratio 6 ' is too small, then tachometer 
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feedback can be used to increase this ratio, if for the specified values 


of zeta and omega, By is negative, etc. 


Now letting H(S) = K 5? in figure (3-8), the characteristic equa- 
a 


tion becomes: 


Se(S) + Pe(S) + KK Ss” + Pk $° +KX¥ S+PRK=0 (3-52) 


By making the same substitutions as in equation (3-49) one obtains: 


k-2 


k k-l 3 
S +...+e, =F e, of + kk Ss 


S + (ae + e-ps +(e 1% te, 


-2) 


2 
(e, or en of +k of )S + (e, A tet es + of (K + e) =Q (3-53) 


It then follows readily that: 


B, = -(K te.) + (e, + KK )w -ewU, te t(-1) Ze why 
a Cla. 
Cc, = 0 
D, = ew - (estKK )w teeet(-1)S ewe U (3-53a) 
+(-1)te wu + (-1) "WG, 
By = penn | (e, + KK,)wU, = egw Ute t(-1) ew 7u + 
yw sy 
Cy = -Kw 
D, wane + ewU, - (e,+KK wu, tet fe wu + 


k-1, kel k 
(-1)Fe wh tu + CD Sw, 


After writing equations (3-53a) in terms of the primed quantities, which 
correspond to the situation of cascade compensation only, it can be seen 
that the following expressions result: 

2 2 


B. = B! + kK w B, = B' + KK wwU 
a 2 a 


1 1 Z 2 
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ny ae on ee = t ut 
Cc, =C, = 0 C, = C, (3-54) 


~ KK Pe D cue AK wu 
1 a a 


a 
H} 
o 
i 
oO 


Solving for alpha and beta one obtains: 
f= -(0! - KK Ww) /(B4 < KK w) (2255) 
B = GB, + KK wl )(D! - KK Fy - (BI + KK w7)(D! - KK art ) 
ae 2° * I a 1 a 2 a 3 
-wKk(B1 + KK _w*) 

There is no straight forward way of observing from equations (3-55) 
how “© ' and 6 ' are modified by acceleration feedback. However, assum- 
ing that one would logically try cascade compensation alone before attempt- 
ing combination, one could easily compute the alpha and beta with accelera- 
tion feedback added since the quantities By» Bo» Di» and D, in equations 
(3-55) would already have been computed. 

Example 3-/7 
Problem: 

Design a cascade compensator with tachometer feedback as shown in 
figure (3-9), so as to place a pair of characteristic roots at zeta =.5 
and omega = 1. K should be 50. 

Solution: 

Referring to example (3-6) it is seen that this is the same problem 
except that tachometer feedback has been added. As was found in example 
(3-6) the value of beta using a single section compensator was too small. 
The problem now is to see if tachometer feedback will increase this value. 

As was determined previously, K = 500, BS = -503, C, = 0, Ds = 9, 


By = -9, Ci = -500, Di = 6, ©' = .0179, and 8' = .0117, where the primes 


2 Zz 


have been added to indicate cascade compensation only. Using equations 
(3-55) one finds: 
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B = (OK, + 9K, + 500K0)/(-503) + K, + .0117 (3-55a) 


It follows from the above equation that for beta to increase the following 
inequality must hold: 
K, > (18K, + K2)/503 
or (3-56) 
Ke < (485K, )/500 
Since one excludes values of K. less than or equal to zero, equation (3-56) 
reduces to: 
0o< K i 
K. can arbitrarily be taken as .4. Beta can then be obtained from equa- 
tion (3-55a) and is found to be .238 which is in the acceptable region of 
a= B75 = 10. 
Calculating alpha from equations (3-54) one finds that: 
Be KK w"/B! ia Se 
The cascade compensator becomes: 


.238 (S + 1.744) , which can be synthesized as an 
(S + .414) 


R-C lag network. The compensated system's characteristic equation is as 


follows: 
st + 8.415S> + 220.387 + 219.068 + 211.69 = 0 (3-57) 


The roots of equation (3-57) are: 


r, = --499 - 5.866 
r= -.499 + 5.866 
r, = -3.708 - 414.077 
r= -3.708 + j14.077 
The roots r, and r, are the desired ones, and since .499 << 3.708, they 


are also dominant. 
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Derivative signal enclosing a cascade compensator, 
In figure (3-10), G = K/e(S) and Gol ¥S + P , where these 


S + P 
quantities have been defined previously. Let H(S) = KS. The forward 


path transfer function is then seen to be: 


gus SS d K ¥ (Ss + P/¥) (3-58) 
ee oe H(S)G_G (S + P)e(S) + K 8 (S + P/x )H(S) 


Since the D. C. gain of G. is unity, the error coefficients for figure 
(3-10) are the same as those given in table (3-1). 
By expanding equation (3-58) into the characteristic equation one 


gets after letting alpha = P and beta = Y¥ 


Se(S) + Pe(S) + K¥ KS” + KPK,S + K 8 S + KP = 0, or after expanding: 


k k-l k-2 2 
S + (AH + ee» + Ce _yA+ e-2)5 +...+(e, + KK. B + Pe, of )s 
+(e, of 1 os KK, of +K Q + es + of (K + e,) = 0 (3-59) 
Then as before one obtains: 

- 2 k-2 k-2 kel k-l 
- 2 
Cy = KK LW 
- 2 k=-2 k-2 kel k-l 
D, = ew +...+(-1) e- U3 + (-1) ei" U9 
k 
+(-1) wy 
s 2 k-2 k-2 k-l k-l 
By = -(e, + KK, ) wt enw Upt...+(-1) e-1" U9 + (-l)° ow Uy 
2 
C, = -Kw + KK LW U, 
D, = -e we +e w*U + Gi ce oar + pre. why 
2 oO 1 2 m-2 k-2 m-l k-1 
+ (-1) Sw, 


With primed quantities corresponding to the case of cascade compensation 
only it follows that: 


— t 
ey 
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l l 
- n! 
Dy Dy 
iret. 2 
ree EC (3-60) 
2 
= $ 
C, = C, + KK Ww U, 
D =D' 
2 2 
and 
of = KK.wpD! - (C! + KK wt,)D! 
C2 Z t" 248 (3-61) 
2 2 
¢ $ 2 oe 
B, Cc, 4 KK,w U,) (B, KK w) KK w 
= to t _ pint 
= (By - KK.w)D) BiD, 
2 2 
i¢nt si ‘ | 
BCC, “ KK .w U,) (B, KK w) KK w 


In figure (3-10) letting H = Key the characteristic equation becomes: 


k k-2 


k-1 3 
S + (Ll + e-5 4. (e.1% 4. e 7-299 +...+(e, +e, + KK, B)s 


3 


2 wd 
+ (e, “+ enh -+ KK Ss” + (e,% + K B + e.)s +o (K + e.) = 0 (3-62) 


Then in terms of the primed quantities it is found that: 


2 
— Rp! 
By By + KK Ww 
= ' - 3 — <= 3 
Cy Cy KK LW U, KK w U, 
= D! 
1 (3-63) 
3 
—- R! 
B, B, + KK w U, 
C, = C! - KK weU 
2 2 a 3 
- 1! 
D, D 
and 
of = -KK wU Dm = (Ci - oR wu )D! 
a ., 202 2 a 3" 4 (3-64) 
(B! + KK wy (Cc! - KK wu >) + (Bi + KK wl )KK wu 
ih a 2 a 3 2 a 2 a 2 


38 





Z 2 
= ne = ' ' 
e (Bs + KK w U,)D; (By + KK ow )D, 
(B! + KK roe - KK wu ) + (B) + KK wall )KK oo 
1 a 2 Big 3 2 ac G2 a2 
Comparing the expressions for alpha and beta for the case of the deriva- 
tive signal not enclosing a cascade compensator to the case of the de- 
rivative signal enclosing a cascade compensator one finds that the latter 
are considerably more complex. So if one had a choice between the two 
methods, the former could be tried first since it is easier to analyze. 


If the values for alpha and beta thus obtained were still not acceptable 


then the latter method could be tried. 


3-2 Dominancy of the specified roots. 

In the preceding examples nothing was done in the calculations to 
make the specified roots a dominant pair. As was mentioned in section (1), 
being able to predict a system's response on the basis of the location of a 
pair of complex roots was based on the assumption that the magnitude of 
the real part of the primary or specified roots was much less than the 
magnitude of the real parts of all the other roots of the characteristic 
equation. In most cases, if the real part of the primary roots is one 
half to one fifth or less of the real parts of all the secondary roots, 
the system is said to be dominant in the primary roots. In many cases 
the system will still meet the specifications even if two pairs of complex 
roots have the same real part, providing the zetas for both pairs of roots 
meet the specifications, and the undamped natural frequencies are such 
that the component time responses are not highly additive. The presence 
of closed loop zeros will also greatly affect the dominancy factor needed. 
For instance even if there exists a characteristic root whose real part is 


closer to the origin than the real part of the primary root, the presence 
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of a closed loop zero could make the residue of the close in root negli- 
gible as compared to the residue of the parimary roots. However, if possi- 
ble, one tries to make the real parts of all secondary roots as large in 
magnitude as possible. 

In the preceding examples it should be noted that in many cases 
there were actually three and sometimes four variable parameters. For 
instance, the forward path gain was usually set at a fixed value in the 
computations so as to meet the minimum steady state accuracy requirements. 
There is, however, usually no reason why the gain cannot be raised above 
the idiom value, thus permitting a third degree of freedom. When cas- 
cade and feedback compensation are employed simultaneously, the forward 


path gain and tachometer gain become the third and fourth parameters. 


3-2-1 A method of employing a third parameter. 
Recall that the system characteristic equation has the following form 
mM 
k 
£(S) = S aS = 0, where a = bot + c. 6 + di. (equations (2-1) and 


K=O 
(2-8)). In order to meet the system specifications, one places a com- 


7 ne | 2 : J’ 
plex root pair at S = - 71") + iw, 1- r 4 , which implies that 


s7 4 2 z w,S + w 0. (3-65) 


Since the coefficients of equation (3-65) are known, this quadratic can 


be divided out of the characteristic equation, leaving a polynomial which 
contains all the secondary roots of the characteristic equation. Since 
only two of the degrees of freedom or variable parameters were used in 
fixing the roots of equation (3-65), the remaining variable parameters 
will appear in the coefficients of the quotient polynomial, and it is 
these coefficients that can be varied to achieve dominance. Instead of 
division to find the quotient polynomial, coefficients of like powers will 


be equated to achieve a set of equations. Let the quotient polynomial be 
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given by: 

£,(S) = i fs = 0 (3-66) 
where n = m-2, i.e., equation (3-66) is order two less than the character- 
istic equation. Using equations (2-1), (3-65), and (3-66) it is seen 
that: 


h 
2 2 ae k 
(s~ + 2 7S + )( Zi f,S) = 2 a, S (3-67) 


a7 = fy.p t+ 2 7 Wy 


2 
ae ine * fan? 71"%4 


: 9 (3-68) 
ay = £ + 2 711 Fy + £1 
2 
ay 2 yee + fe) 
im 2 
Bo r or 


The formulas (3-67) can be solved for the coefficients f in terms of the 


coefficients a. The solution will be made for the following cases; 


Case of k = 3, n=1 


Equation (3-67) becomes: 


2 2 3 2 
(S~ + 2 ame 7) CEs ee”? S + a8 + aS ta, 


Equating coefficients of like powers one obtains: 


eS. = Le f 


~ 
Ki 


pare” © 7 aad 


2 | 
£1 + 2 7 tac 


os 
i 


1 1 


Solving for the coefficients f results in: 


41 





Dee 2 = 
En oie Cl Pee 7) ~ =e (3-69) 


Case of k = 4, n = 2 


Proceeding as before the a coefficients become: 


a, =le=f 
a, = 2 714 + f 


+2 2415, + £0 


my (3-70) 
a, =f — + 2 wif 
i 11 71 1 
_ 2 
a * EY) 
When solved for the coefficients £, equations (3-70) yield: 
f, = ] 
7 7 ' , 2 
f. =a, - 2 7") = a,/2 71" w,/2 ¥ 1 a/2 3 YY 
= 2 
f - 1/w, (a, - 2 ae (3-71) 


rh 
il 


oP . 2 
9 Ao/Wy = ay/2 2 yw, - agwi/2 z,) + wy) = a, - 2 7 ees 


2 22 
w) +4 Vee ees + wu, 


Case of k = 5, n = 3 


Proceeding as before the coefficients a are: 


a =le=f 


5 3 
ee Fa) tube 
_ “ape 
a, = WwW, BB? Paine +f, 
a, =f +2 z,w,f, + 2 
2 an i 
a, = 2 7 123', f +f ee 
= Z 
=o 7 cael 
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Solving for the coefficients f one obtains: 


f£, = 1 


f= lay) = 2 Fa, /y, 


c sa, - 2 3,0, = a,/wy = ais 74 wW 
f = a,/2 7M, - a,/2 3, + awe -w/2Z, 


2 
1 


+s 3 yw, 


£ awe - 23 4a fw, =a, - 24 »¥) 4, + Uw (3-73) 


3 
f =a hae =a_/2 w, - a.w,/2 ea 
OMe BSlesce a | Soe aan) 
leas ies 2 4 Wot 
In formulas (3-69), (3-71), and (3-73), the coefficients a are of the 
form a, = by of + Ch B + d,- In section (3-1), formulas for alpha and 
beta were developed for various forms of compensation techniques. Since 
the formulas for alpha and beta will in general contain other variable 
parameters, then the coefficients f as derived above will be functions of 
these parameters. In most cases the coefficients f will be functions of 
only one parameter (or they can be made to be). Hence the roots of £,(S) 
can readily be placed algebraically for the cases of n= 1, and n = 2. 
For the case of n = 3, a root locus sketch will usually suffice. Although 
the coefficients f have been derived for only up to a fifth order case, 


they could easily be obtained for higher order cases if necessary. 


3-2-2 Applications of the dominancy technique. 
Example 3-8 (Third order characteristic equation) 
Problem: 
Compensate the system of figure (3-11) with a cascade compensator to 
obtain: 
1. Characteristic roots at zeta = .5 and omega = 40. 


De K = 250. 
e 
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3. The specified roots are to be made dominant. 

Solution: 

The characteristic equation of figure (3-11) is: 

6G Pape + ore R © )S + Keone 

or s3 4 (4 ol 8 + (4L +KB)S + Kx = 0 

where alpha = P and beta =Y¥ 

Now G = K/e(S) = K/(S* + 48), soe = 0, e,=4,e,=1,U, = 1, U, = 0. 

Equations (3-42) are now applied to obtain: 


B. = -—K+w = -K + 1600 


1 
C, = 0 
D, = ass? ~ we = Bean 
1 2 
. 27 
B, = -4Awt+w U, = 1440 
Co = -wK = -40K 
= hit J 
D, = 4wU, - wU, = 6400 


From equation (3-43) are obtained: 


oh = 5.76x10" 
-K+1600 
6 = 1440) -5.76x10°) - (-K + 1600 6400) (3-74) 


-4OK(-K + 1600) 
Now for any value of K, equations (3-74) will provide a value of alpha 
and beta to provide characteristic roots at zeta = .5 and omega = 40. 
The value of K will now be chosen on the basis of dominance and steady 
state exten considerations. To satisfy the error specifications it is 
necessary that K 2 1000. Since £,(S) is of order one, the appropriate 


equations are (3-69), hence: 


A 





TI-€ oainetA 
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ia 2 Foil 
0 7 Af) = (ay > wy)/2 3 wy) = ay 29 144 


hh 
ti 


From the characteristic equation it is seen that 


a, = 4 +0 
=4o0 +K 
a, B 
a = Kol 

Oo 


The real part of the specified roots is 71"1 = 20. Arbitrarily choos- 
ing a dominance factor of five, the dominance criteria becomes: 


fo > 5 g i") = 100. To satisfy this requirement, the simplest form of 


EO will be chosen, namely EO = al lw, Therefore it is seen that = 
K o£ /1600 = (5.76x10") = 36K > 100, wh@eevemidernmce) 


1600(-K + 1600) (1600 - K) 
was employed. 


The above inequality implies that K 7 1180. Since K 71180 also 
satisfies the error specification, a value of K = 1200 is chosen arbi- 
trarily. Using this value of k, the following quantities are computed: 
alpha = 144, beta = 4.2, and EO = 108. As a check, the expression EA = a, 
-2 3 Wy can be employed. Hence: 

tO = (4 + 144) - 2¢(.5)(40) = 108. 
Example 3-9 (Fourth order characteristic equation) 
Problem: 

Compensate the system shown in figure (3-12) employing tachometer 
plus acceleration feedback to obtain: 

1. Characteristic roots at zeta = .5 and omega = 2. 
2 Woese2. 


e 


3. The specified roots should be dominant. 
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Solucion: 
The characteristic equation from figure (3-12) is; 


3 
s+ 4+ 16.58 + (73 +o )S* + (82.5 + B )S425+K = 0 


Here alpha = RK. and beta = KK 


3 


G = K/e(S) = K/(S* + 16.5S> + 738° + 82.58 + 25) 


By inspection it is seen that: 


aa 25’, e, = S255 e, = 73. e, = L635; e, = ie U, = ], U, = 0, 
and U, = -l,. By employing equations (3-11) one can find: 
alpha = (185 + K)/4 
and 
beta = (K - 23)/2 
From the characteristic equation it is seen that: 
a, = 1, a, = 16.5, Bp Pee a F B2 eo pe eet K. 
Since the quotient equation £, (Ss) is a quadratic, i.e., 5? te f,5 + 


fo = 0, it would be desirable from the dominancy standpoint if {7 5 7 ea 
= 5, However, looking at the dominancy equations for this case (equations 


(3-71)), it is seen that one of the several expressions for ff is f aos 


2 71: Since all the expressions for f, have to be simultaneously 


i 


satisfied, it is seen in this problem that f, is a fixed constant since 


a. and 714 are fixed. This is true even though the remaining express- 
ions for fy involve one or more variable coefficients. Therefore fy is 

found to be 14.5. Observing the pertinant expressions for f it is seen 
that none of them contain only constants, so the simpler expression , = 


1 1 


exists, but the system's performance can be further improved by choosing 


2 
a /w, is chosen. Now since f, = 14.5 > 5, a dominant situation already 
O 


a reasonable value of zeta and omega for the secondary roots. From the 


error specification it is necessary that K/25 = 12, or K 2300. Now 
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£,(8) = So ieee awe or £,(S) = S* 4 14,564 6525 Connon 

K = 300, £ (S) becomes: 5? + 14.58 + 81.25. Therefore, 2 a o> = 
l 

aj = 81.25 or Wy = 9. Then zeta = .806. These are reasonable 


values for E 9 and Wo since the secondary roots taken by themselves would 


14.5, w 


produce much less overshoot and a much smaller settling time than the 
primary roots. Using this smaller value of K one can compute alpha and beta. 
alpha = 121.2 beta = 138.5 
Since alpha = KK, + 300K, and beta = KK, + 300K. then K, = .405 and 
K = .462, 
a 
As an added bonus of the method, all roots of the characteristic equa- 


tion are now known and the time response could be computed if desired. 


Example 3-10 (Fifth order characteristic equation) 
Problem: 
Compensate the system shown in figure (3-13) to obtain: 
1. Characteristic roots at zeta = .5, and omega = 4. 


2. kee S 
e 


3. The specified roots should be made dominant. 
Solution: 


The characteristic equation is after letting alpha = P and beta = ¥ 


Pe GR las Glo 4ueaies = (14a: KK, + S4oL)s* + (148h + 


KK, of + K 3+ 80)S + of (K + 80) = 0. 


3 2 


G = K/e(S) = K/(S* + 178° + 848° + 1488 + 80) 


From the. above it follows that: e, = 80, e,= 148, 2 = 84, e, = 17, 


= 0, U, = -l, and U. = -1. Employing equations (3-42) 


e, = 1, U. = 1, U 4 5 


4 Z 


one obtains: 


3 
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B, = K + 1196 B. = 492 


1 2 
Cece Co eas 
D, = -Uae6 D, = -1276 
Therefore: 
alpha = 1986/(K + 1196) beta = (1276K - 555000)/[4K(K + 1196)] 


For beta to be positive it is necessary that K be greater than 435. The 
ee cacy specification implies that K be greater than 640. For K = 640, 
beta = .0555 which is out of the desired range of .l1 = 6 S108 
Of course, increasing K makes beta even smaller. At this point one could 
design a multiple section compensator in accordance with section (4-6-2) 
or feedback compensation could be added. The latter course of action 
is chosen and tachometer feedback not enclosing the cascade compensator 
is chosen. The applicable equations are (3-54). From them are obtained 


the following, where the above alpha and beta now become “' and 6 '. 


_ 1 (1986 - 16KK.) 
alpha K +1196 ‘ 
‘ 2 
beta = -492(4)K, - 1986K, + 16KKE oK, +B 
K + 1196 


To simplify the analysis let K = 640 which meets the accuracy specifica- 


tion. Then alpha and beta become 


alpha = 1986 ~ 10210K, aii 
1836 
Here alpha is positive for K. less than .193. 
Geta = 5. 64K. + 1.83K, + .0555 (3-76) 


The appropriate dominancy equations are (3-73) where it is seen that 


neither f. f., or f, are restricted to constant values. From the charact- 


1 2 


eristic equation it is seen that: 
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a, = ole 
a, = 17 + 84 

= of 
a, 148 + KK, + 84 


1480L + KK d+ K 8 + 80 


jah) 
ii 


oh (K + 80) 


ja) 
it 


On the basis of these equations the following expressions are chosen 


from equations (3-73) for the £, (S) coefficients: 


= - = of 

f. ay, 2 y "1 +135 
= 2 = 

5 = a., 2 y 1134 + U,w = 130% + 16 
ll 2 ee 

f, = ai/w, = 1804 


Then £,(S) becomes: 


s> + (of + 13)8* + (134 4+ 16)8 + 180% = 0 (3-77) 


Equation (3-77) contains only one variable so it can be put in root locus 


form as follows: 


s> + 138% + 168 +L(87 + 13S + 180) = 0 

After dividing by 53 oP 138° + 16S the above equation becomes: 
(s* + 138 + 180) = (S + 6.5 F jll.7) = -l (3-78) 
S(s* + 13S + 16) S(S + 1.37)(S + 11.6) 


By plotting the root locus poles and zeros of equation (3-78) or 
by inspection if one wishes, it is apparent that there will be a real 
root of £,(S). Also the magnitude of the real part of the complex roots 
can never be greater than 6.5, and it approaches this value as alpha tends 
to plus infinity. For dominancy therefore, a large value of alpha is de- 


sired. From equation (3-75), the maximum value of alpha is for n= oe 
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but then beta is too small since Ke = () implies the case of cascade conm- 
pensation only. It is obvious that K. should be as small as possible, 
which from equation (3-76) implies that beta should be as small as possi- 
ble. Applying the lower limit of beta = .1 to equation (3-76) one can 
show after solving the quadratic in KL that K. = .023 is the necessary 
value. 

From equations (3-75) and (3-76) one can solve for alpha and beta 
obtaining: 

alpha = .955 = P beta = .1 = ¥ 
The cascade compensator is then: 


1(S + 9.55) 
(S + .955) 


This can be realized by an R-C lag network. Equation (3-77) for £,(S) 
becomes: 
s? + 13.9558" + 28.45 + 172 = 0 (3-79) 

A root locus of equation (3-77) for 0 = “<q is shown in figure (3- 
14). The root locus was done by the computer program presented in sec- 
tion (5-2). The roots of equation (3-79) were obtained from the printed 
out data from the program and are as follows: for alpha = .955; S=- 
-585 + j3.64 and -12.8. Since .585 is less than 71 1 Which is two, the 
primary or specified roots are dominant. Now the maximum value of alpha 
is obtainable when K = KR = 0 and even for this extreme case, the speci- 
fied roots will not be dominant. To solve the problem, therefore, a 
compromise on the specifications must be made. One approach would be to 
specify omega as some value less than four and repeat all the preceding 
calculations, but one could also do the following: 

From equation (3-75) observe that alpha can be decreased by in- 


creasing Ke: From the root locus data it is observed that if alpha = 


a3 
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117, then ? 5 = .5, and Wo = 1.34. If this reduced value of omega 
can be tolerated then the problem is solved since the secondary roots 
located at 7 9 and W, are now dominant and only the settling time has been 
modified (due to the decreased value of omega). On the other hand since 
the roots 71 and w, are still located at zeta = .5 and omega = 4, the 
overall system response will probably be acceptable to the designer. 
The final results are therefore as follows: . 

K, = .173 obtained from equation (3-75) 

alpha = P = ,117 

beta = gamma = .5405 obtained from equation (3-76) 
Characteristic roots are located at zeta = .5, omega = 4; zeta = .5, omega 
= 1.34; and a real root at S = -11.8. The cascade compensator is of the 


forms: 


-9405(S + .216) 


(S + .117) 


This can be synthesized as an R-C lag network. 
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3-3 Some sketching techniques. \ 

The graphical solution will be discussed in following sections and 
of course this involves plotting curves. Due to the complexity of the 
parameter plane equations, two methods of curve plotting are recommended. 
The easier and faster method is by means of a digital computer. A com- 
puter program which will do this is presented in section (6). The second 
method is by means of a desk calculator or slide rule. To facilitate the 
use of the latter method, the parameter plane equations are manipulated to 
put them into forms more suitable to non-computer plotting. 

Due to the time limitations usually imposed on digital computer 
operation and due to the time consumed in plotting the curves by hand, it 
would be desirable to just sketch, say the zeta equals zero and the zeta 
equals one-half curve for various values of omega to see if the type of 
compensation proposed will do the job. This type of rapid sketching is 
also helpful when choosing a scale for the digital computer graph plot. 
For this reason, in the following sections, special characteristics of the 
@ta equals zero and the zeta equals one-half curve have been derived for 


characteristic equations of order two through five. 


3-3-1 Table of symbols. 

In table (3-2) are listed various symbols which apply to the deriva- 
tions made in section (3-3-2). The change of variables as indicated in 
table (3-2) is useful because many of these products and sums appear re- 
peatedly. When doing hand computations of the parameter plane curves it 
is only necessary to compute these quantities once, so substantial labor 


is saved. 
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“boc, + bic, 

b cy - b,¢ 

cd, - Cod, -<¢ 
- - c,d, 

b,¢ = bod, 
bod, ~ b,d, + b, 
dco - doc, 
boc, - bofo 

b d, - bod. 
boc, - b,c, 
bac, - b C4 
Cod, = c,d, -c 
Cad | - Cc d., 
bd, = bod, +b 
bod, = b,d, 
B+ J 

D+ Lt CO 

P - bo + F 


d. -c.d.~.+c 
oO 


Som 1°3 


TABLE (3-2) 
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ul 1} il Nl 


ner oe © 80 2 SO Oo FF HPO 
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é 
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3-3-2 Basic derivations. 
Case I (Second order characteristic equation) 


The characteristic equation is of the form: 


2 
¢ X -1 C of = 
Sv + (b,%+ c,6+ d,)s Dot cbt d 0 


Employing equations (2-10) one can obtain: 


By = oe B. = “bow 
C, = aS C, = -c,* 
id tie Dis ede & Gal 
1 oO 2 1 2 
Solving for alpha and beta using equations (2-11) results in: 
alpha = (cw - cola + M)/K (3-80) 
2 
beta = (-b,w + b Uow + N)/K 


To find the maximum and minimum points, the first derivatives of alpha 


and beta with respect to omega are taken and set equal to zero. 


d & /dw 


(2c,w = cU,)/K =Oorw= c U,/2c, (3-81) 


d 6 /dw = (-2b,w + b.U,)/K = 0 or w = bU,/2b, 
In equations (3-80) after letting w = 0 one obtains: 

of = M/K 6 = N/K (3-82) 
Letting w tend to plus infinity one obtains: 


alpha —> +00 beta —> + (3-82)° 


Equations (3-80) through (3-82) are valid for all values of zeta 


between zero and one. Since U, = 0 when zeta = 0, then equations (3-80) 
and (3-81) become: 
agree = (cw + M)/K beta = (-b,w + N)/K (3-80a) 
d“/dw = 2c, w/K w = 0 (3-81a) 
dB /dw = -2b,w/K w = 0 
3 


The relative magnitudes of the coefficients determine whether alpha 
and beta approach plus or minus infinity. 
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Equations (3-82) remain unchanged. Also since U, = 1 when zeta = .5 


one can obtain: 


aan (c,w -c.w+M)/K beta = (-byw" + bow + N)/K (3-80b) 
d & /dw = (2c,w - co) /K w= c /2c, (3-81b) 
d © /dw = (-2b,w +b )/K w = b /2b, 


Equations (3-82) remain unchanged. 


Case II. (Third order characteristic equation) 
The characteristic equation is of the form: 


3 n 
So + (b,L+ Cy B+ d.)s + (b, ol + c, 6+ d,s + bot + c b+ d. = 0 


proceeding as before the following expressions are obtained: 


4 2 3 2 
Ripka = W (-c,U, + c5U,) ~ c,Uow + (L+ Cy + cU,)w + U,Rw +M 
ame 7 U,tw +K 
(ab v2 + b,U.) + b,U = + (P-b -bU Neh + U,Uw + N 
beta = 2.2 2a 1 2 Oo o 3 2 
Tar + U,Iw +K 

For w = 0, 

alpha = M/K beta = N/K (Same as equation (3-82)) 
For w>+ 0% 

alpha = +0 beta = +0 (Same as equation (3-82)) 


Due to the increased complexity of the expressions for alpha and beta, 
the derivatives are only computed for the zeta equals zero and the zeta 
equals one-half curves. 

Letting zeta = 0, equations (3-83) become: 


4 2 4 2 

niptoc oe YM ae= oe a (3-83a) 
jee + K St ih 

For w= 0 and w= , equations (3-82) apply. 


Taking derivatives of equations (3-83a) one obtains: 
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4 2 
a? of Vein if W CJ +w 2c ‘ + KL - JM 


2) ees ee Oey eee ee 


Gar + K)? 


INO 


The derivative goes to zero for: 


2 
we . ~2c,K a (2c,K) ~ 4c5J(KL - JM) 


2c,J 


4 2 
a2 B few? = W(-bjJ) - w'(2b,K) + KP - NJ 


(Sw" + K) 


The derivative goes to zero for: 


| 2 
wee 2boK + (2b,K)" + 4b,J(KP - NJ) 


2boJ 


Letting zeta = .5, equations (3-83) become: 


Z 
alpha = caw" = ew + (L + cw + Rw + M 


in + Iw + kK 


beta = -byw" + bw +(P - bw + Uw+N 
2 


Iw +Twt+kK 
For w = 0 and w =o , equations (3-82) apply. Taking derivatives 


equations (3-83b) one obtains: 


4 
T + 2c, J)w + (4c 4K = 


5 
fh = J - 
d /dw [2c Jw + ( 3c, J + 4c 1 Co 


2e,T)w + (JR + LT +e T = 3¢,K ~ 2IR)w" + (-2JM + 


2LK + 2c K) + KR - MT] /(Jw? + Tw + K)? 
4 5 4 3 
dh /dw.= [-2b,Jwo + (-2b,T + b,J)wo + (-4b,K + 2b,T)wo + 


(3b,K + PL - bT - UJ)w? + (2PK - 2b K + 2JN)w 


+ (UK - NT)]/(3w* + Tw + RK)? 


(3-84) 


(3-84) 


(3-84) 


(3-84) 


(3-83b) 


(3-83b) 


of 


(3-85) 


In the general case, equations (3-85) involve fifth order poly- 


nomials in omega. In specific problems, however, some quantities 


in 








these equations will be zero, enabling one to more readily solve for the 


critical values of omega. 


Case III. (Fourth order characteristic equation) 
The characteristic equation is of the form: 


4 3 2 
S + (b, A +c, 6 +b d,)S + (boa + Cy 8 + d,)s + 


(be +e, PB td)Stb ete Atd, = 0 
Proceeding as before the following expressions are obtained: 
a 6 2 5 4 
alpha = [ w c,(U, - UU,) +w cy (U,, = U,U,) +w (-C - cy 


2 3 
-- c,U, + Us (C + c,)) +w (U,& - co) - c U,) + 


2 
Ww (-U,D +Lt+ Cc.) + wU,R +MJ/A 


"2a 2 5 (eee 
beta = [w b, (UU, - U3) +w b,(U,U, = U,) +w (UA (E - 
3 
by) © ty. UE b,)) +w u,(Z + Boe + b U,) 
2 

‘vet D - U,F) + wU,U + N/A 

ttt. 3 2 

A=wA(U, - U,) +w YU, +w (-BU, + J) +wIU, + K 


For w= 0 andw=© , equations (3-82) apply. Letting zeta = 0, 
equations (3-86) become: 


ame? ecg Gee 





alpha = S 
ne + eur + K (3-86a) 
6 4, 2 
Pun & -b.w + Ew +iIw +N 





ae + Gw" + kK 


For w= 0 and w = ©, equations (3-82) apply. Taking derivatives of 


equations (3-86a) one obtains: 


a of Teak = [w Ac, + w°3e 46 + w (AH + CG t+ 3Kce,) + ww (2CK + 


2 


GH - 2A(H + M)) + KH - G(H 4m) ]/(Awt + Gw* + x)? 


(3-87) 
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ne B eae = [-w°Ab, - w'2b,6 + w' (EG - AI + 3Kb,) + a (GI + 


2 


2QEK - 2AN) + KI - NG]/(Aw + Gw* + x)? 


In the general case, equations (3-87) involve eighth order poly- 
nomials in omega. In specific problems, however, some quantities in 
these equations will be zero, enabling one to more readily solve for the 
critical values of omega. 


Letting zeta = .5, equations (3-86) then become: 


6 5 4 3 2 
alpha = °3” ~ So¥ * (C + c,)w + Xwo + (L+c wo + Rw + M 
pace + yw? + aa + TwtkK 


(3-86b) 
6 5 4 3 2 
ae -b,w F bow + (E - b,)w + Zw + (P - b)w +UW+N 


ee + Yw> + jee + Iw+kK 
For w = 0 and w =© , equations (3-82) apply. Due to the increased 
complexity of equations (3-86b), it does not appear practical to compute 
the derivatives for the general case. 
Case IV. (Fifth order characteristic equation). 
The characteristic equation is of the form: 


5 4 3 
Ss" + (b, of + Cy, B + d,)s + (b, 06 “F C4 B+ d,)s + (b, of + cy B 
Z 
+ d,)s + (b, o +e, @+4))s aon te B+d, 0 
Proceeding as before the following expressions are obtained: 
alpha = mie (UL + u*) - ne (U.U. + U,U,) + wi (U U 
4° °5 4 2 2") 3°4 24 
2 5 
=U c, (U, + U,U,)) +w (U,U, K + Ue - c,U,) 
+0 (U2 m- U,m+ UW + c_) tw O+u A) + 
2 3 3 Oo 4 or 


w (U4 7+ L +c) + WR + MI/A (3-88) 


2 
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ee ee 





8 2 7 | 6 
beta = [-w b, (U, - UU.) + Ww b,(U3U,, + UU.) - w C (ULU, 
yz - b.(U.U, + U.)) + w((UU, - U,)S + b,U,) + 
a as pan oh WY, 
wH(-U.(E = b.) +U,T + UE - by) +b.) tw 
3 1 3 2 1 o 5 


Z 
(U, A+ U,( + b)) + ow (-U,F +P - b)) - wu 5U +N2/A 


® 4G 2 5 4 2 
Aw, - Uz) +H 6 (U, - ULV.) + w ¥ (UU, - UZ) + 


3 2 
W (U,€ + U.¥) +wi(J- U,B) + wu oT +K 


For w= 0 andw=© , equations (3-82) apply. Letting zeta = O, equa- 


tions (3-88) become: 


8 6 by 9 


41 (3-88a) 


8 6 4 2 
beta wb, -w P (bo +1) +w (EB - dy —P +b) @ Ria beaee 


41 


4D, ais - oer ot Tea) +K 


Letting zeta = .5, equations (3-88) become: 


alpha = [wc = we ( + + 2c.) + wc, 5 ce. ) + wth 7 wR : 


B)+wWtc.) + wR + Ml/ 42 (3-88b) 
beta = [-w'b, ~ w'b, + e (2b, - 1) + w (S-b,) ae 
3 Z 
(E - b, -b) +w (Z +b. -MK)tw(P -bi) +wtn/A, 
mn a 6 -wo-wV¥+w(-6€)+wJt+w+K 
il | 
Example 3-11 (Case I example) 
Problem: 
Sketch the zeta = 0 and the zeta = .5 curves for the following 
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characteristic equations. Let the abscissa variable be alpha and the 
ordinate variable be beta. 

s7 + (f+ 2)8+B4+1=0 
Solution: 


From equations (3-80) it is found that: 
2 


alpha = Upw -~ 2 beta =w - 1 
From equations (3-81): 
SoU d 8 /dw = 2w 


From equations (3-82): 

alpha = -2 and beta = -1l for omega = QO. 
Therefore for zeta = 0 the following relations apply: 

alpha =-2 beta = ne - 1 and d 6 /dw = 2w, 
which implies a minimum for beta at omega = 0. The zeta = O curve is 
therefore a vertical line in the alpha-beta plane, going from alpha = -2, 
beta = -1 to plus infinity as omega increases from zero to infinity. 


Setting zeta to .5, one finds that: 


- | 


it 
<= 


alpha = w - 2 beta 


or 


walk +2 beta = (ol+ 2)? - 1. 
Tne alpha-beta curve is therefore a parabola with vertex at alpha = -2 
and beta = 0; symetric about the line alpha = -2, and opening in the plus 
beta direction. The curve can readily be plotted after calculating a 
few critical points. 

The above curves are plotted in figure (3-15). 
Example 3-12 (Case II example) 
Problem: 


Repeat example (3-11) for the following equation: 


s> + (kt 187 + ( B + 2)8 +3 = 0 
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Solution: 

From equations (3-83a) for zeta = 0, it is found that: 

alpha = -1 + an" beta = ee - 2 
For omega = 0, alpha = plus infinity and beta = -2. 
From equations (3-84): 

a> wees ae a7B Jaw? =4 
But d 6 /dw = 2w, so alpha has a minimum at omega equals infinity. Beta 
has a minimum at omega = 0 and it is -2. The curve is asymptotic to the 
lines alpha = -1 and beta = -2. 


The following points are readily obtained: 


omega alpha beta 
1.414 a 0 
1.782 0 1 


From equations (3-83b) for zeta = .5, it is found that: 


alpha =w-1+t 3/0" beta = — -2 + 3/w 


For omega = 0, alpha = +© , and beta = +00 
Also, 


d © /dw=1 - an d B /dw = 2w - 3 fw" 


It is therefore obvious that: 

ad B/do = (20> = sige - 6), from which the critical points are 
found to be: w = 0 and w = (3/2) 1/3 = 1.145. This point is a minimum 
Since the second derivative is negative. 


The following values can be obtained: 


omega alpha beta 

1.145 2.435 1.93 
ad 300 28 
1 3 2 
3 2253 8 


The curves are sketched in figure (3-16). 


Example 3-13 (Case III example) 


= Sie) ei) § tee Oe 








Problem: 
Repeat example (3-11) for the following equation: 


5 + (eK + 1)s> + 25° 


+S 6 eee 0 
Solution: 
From equations (3-86a) for zeta-0: 
alpha = eae - 1) beta = (aa + aw + peayy cee - 1) 
As omega Pend: to plus infinity, both alpha and beta tend to minus 
infinity. As omega tends to zero, both alpha and beta tend to zero. 
Also: 
doh faw? = [-(w? - 1) + ww" - 1)? 
d°B faw? = (Cw? - 1)(-3w! + 6w2 -2) - (-w® + aut - aw/@w? - 1)? 
It can be seen from the above derivatives that d & /d& = 0, implies six 
critical points. It appears unpractical to compute them. 
From equations (3-86b) for zeta = .5, it can be shown that: 
alpha = se + 2w = w(2 - Ww") 
beta = Be ~ oat ~ oe + ow? = er - ow? - w+ 2) 
As omega tends to plus infinity, alpha tends to minus infinity and beta 
tends to plus infinity. As omega tends to zero, alpha and beta both tend 
to zero. 
Using a slide rule or desk calculator along with the above informa- 
tion the curves can be plotted as shown in figure (3-17). 
Example 3-14 (Case IV example) 
Problem: 
Repeat example (3-11) for the following equation: 
5° +38' + 5° + 257 + A~AS+f =0 
Solution: - 
From equations (3-88a) for zeta=0: 


alpha = w? = a beta = ow" - ay 
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exists at omega = .578. The following points can readily be calculated: 


omega alpha beta 
.378 yee2 . 333 
.816 pees 0 
1 0 -1 


From equations (3-88b) it is seen that: 


alpha = -3w? + 2w beta = w - =) + owe 
Therefore: 
ol 7 2 4 2 
d “/dw = -9w + 2 d B /dw = 5w” - 3w + 4Gw 


Hence d 6 /do& = 0, implies that oo - 3w+4) =0. It is not 
worthwhile to factor the above cubic in omega to obtain the three corres- 
ponding critical points. The critical point at omega equal zero indicates 
that the curve starts at the origin. A slide rule or desk calculator can 
now be employed Scag with the above information to plot the curves as 


shown in figure (3-18). 
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Example 3-15 
Problem: 
For the system shown in figure (3-19) let alpha = KK, and beta 


= Ky Ko» and satisfy the following requirements: 

1. Assuming that a computer will plot the curves for zeta = 
0, .5, and 1, determine an appropriate graph scale. 

2. Only first quadrant values of alpha and beta are of 
interest and the graph is eight inches wide and fourteen 
inches high. Alpha is the abscissa and beta is the 
ordinate. 

3. Due to bandwidth considerations, omega should be less than 
1500. 

Solution: 


From figure (3-19) the characteristic equation is determined to be: 


s? + 11008* + (10° +4 )8s +8 =0 


From this it is seen from table (3-2) that: 


L=l P =-1100 J=0 
M = 10° N = 0 K = -l 
R = -1100 T =0 U=0 


Using equations (3-83a): 


alpha = i + M) / (Sw? + K) = a = 10° 


(pwr + N)/(Sw* + K) = 1100w" 


beta 
For omega = 1500: 
alpha = 21.5x 10° beta = 24.7x 10° 


From equations (3-84) or by taking derivatives directly: 


d? &. faw> = 1 ae B fw = 1100 
Hence: 

d 6 /do& = 1100 
Letting zeta = .5 and employing equations (3-83b): 
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alpha = (Rw + 107) / (Tw - Ll) 


beta = fan - 1100w” + Uw)/(Tw - 1) 


Therefore: 
5 2 
alpha = 1100w - 10 beta = -w (w - 1100) 
It can be seen that the zeta = .5 curve is not in the first quadrant for 


omega greater than 550, so this curve does not influence the scaling pro- 
blem. 

The origin point of the curves for omega = 0, is obtained from equa- 
tions (3-82) resulting in: 

alpha = -10° beta = 0 

From the above data it is concluded that it would be best to scale 
using the values of alpha and beta for the zeta = 0 curve corresponding 
to w = 1500. 

Therefore: 

alpha-scale = 3x10" units per inch 

beta-scale = 2x10° units per inch 
Example 3-16 
Problem: 

For the system shown in figure (3-20) it is desired to have roots 
with zeta greater than .5 and a maximum error coefficient. Plot para- 
meter plane curves to find the best values for K and Z. 

Solution: 
The characteristic equation is: 
4 3 2 


S' + 15S” + 150S° + (100 + K)S + KZ = 0 


Let alpha = K and beta = KZ to obtain: 


st 4 15s? + 1508" + (100 + of )s+B =0 
From equations (3-86a) corresponding to zeta = 0: 
alpha = je + Hae + M)/(Aw! + Gow" + K) 


2 4 2 


beta = (me + Iw +N)/(Aw + Gw + XK) 
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Using table (3-2) one can obtain: 


Cla K = -l R =-150 
H = -15 E= 1 Z= 15 
M = 1000 I = -150 P = -150 
A= 0 N = 0 U = 0 
G=0 X = 1 
Then: 
alpha = eae - 1000 beta = ates ~ 150) 
Hence: 


Beta is greater than zero for a less than 150 and alpha is greater 
than zero for = greater than 66.5, so both alpha and beta are greater 
than zero for 66.5 < w~ < 150. 

Letting zeta = .5 and employing equations (3-86b) it is found that: 

alpha = er ~ 15) - 1000 beta = eee ~ 10) 

From these it is seen that for omega greater than ten, beta is less than 
zero and alpha is less than zero. For omega greater than zero and less 
than ten, beta is greater than zero and alpha is less than zero. Hence 
the zeta = .5 curve is not in the first quadrant. Since the order of 
the characteristic equation is only four, it is safe to assume from the 
above data that all constant zeta curves for zeta greater than or equal 
to .5 are not in the first quadrant. Since the specifications of the 
problem cannot therefore be met, it is unnecessary to plot the curves. 
Example 3-17 

Problem: 

Referring to figure (3-21), 

1, Determine a scale for alpha and beta and plot the curves. 

2. Place a pair of roots with zeta = .5. 

3. Make the roots dominant. 

4. Due to bandwidth considerations, omega must not be greater 
than 4000. 
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Solution: 


From figure (3-21) the characteristic equation is found to be: 


s* + 10058° + (5100 + & )s* + (107 +50 )S + 100% 4+ Caan 


Here alpha = KOK, and beta = KK: 
Since the purpose of plotting the curves is to determine alpha and 
beta to meet specifications 2,3, and 4, the sketching techniques are 


first employed to determine if curve plotting is necessary. 


From table (3-2) one can find: 


Coe 0 L = -l R = -5100 
X= 1 M = 10° E = -1000 
> 7 
Z = 4925 U = 5.1x10 N = -10 
A=0 Y=0 J=0 
T=1 K = -5 P = 74600 


From the above values along with equations (3-86b) one finds: 
alpha = i (w - 5100) + 0; Wigs - 5) 

It can be seen that for omega greater than 5100, alpha is positive, 

and for omega less than 5100, alpha is negative. 
peta fv (w = 1005) + 4oau> 4 74500w" 410, 0¢. 051m may 
(w-5) 

Hence for omega greater than or equal to 5100, beta is positive. 

To avoid*positive feedback, it is concluded that zeta = .5 can be 


obtained only for values of omega greater than 5100, which violates 


specification 4. It is therefore not necessary to plot the curves. 
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3-4 Graphical solutions on the parameter plane. 
3-4-1 Advantages of the graphical solution. 

The previously discussed algebraic solutions have the disadvantage 
that a fixed value of zeta and omega must first be chosen. In some cases, 
the remainder polynomial can then be modified to ensure that the specified 
roots are dominant. However, it is not always possible to guarantee that 
roots placed at a specified location can be made dominant, and a trial and 
error procedure may have to be employed to achieve the best values for the 
various parameters. Trial and error may also have to be used in the de- 
sign of cascade compensators where a specified root location may require 
parameter values that are not physically realizable. 

In eee instances, the calculations have to be redone in terms of 
slightly modified specifications or a different means of compensation may 
have to be used. In general, system specifications are not rigidly fixed, 
but can be met by a given range of values or by some upper or lower limit. 

To avoid this trial and error analytical procedure one can employ the 
graphical solution. If a net of curves is plotted by a computer or other- 
wise, one can, by picking an M point or operating point in the parameter 
plane, read from the curves the n roots of the Ace order characteristic 
equation. The trial and error procedure can then be done visually to pick 


an operating point which best meets the given specifications. 


3-4-2 Some examples of the graphical solution. 

In this section only first quadrant values of alpha and beta are 
assumed to be of interest. The graphs. were plotted with the aid of the 
computer program presented in section (6-1). At the end of each curve on 
the graph is a letter and a number. The letters are abbreviations for the 
following quantities: 

Z(zeta), S(sigma or the real part of the complex variable S), ZW 


(zeta-omega or the real part of the complex roots), w(omega or the undamped 
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natural frequency of the complex roots). 


Example 3-18 


Problem: 

For the system shown in figure (3-22) set K at the stability limit. 
Place a dominant root pair within the following region of the S-plane: 
45 %7 2.7, and2 £ w <€ 6. Both tachometer and acceleration feed- 
back can be eed If possible use only one or the other. 

Solution: 

From figure (3-22) the uncompensated system's loop transfer function 
is: 

GH = -1 = K/[S(S + 10) (s? + 5S + 100)], which when expanded be- 
comes: 
5 + 15s° + 150s" + 1000S + K = 0 


To determine the value of K at the stability limit the Routh array is 


formed: 
1 150 K 
15 1000 0 
1250 15K 0 
1.25x10° - 225K 0 0) 
15K 0 0) 


From the Routh array the stability limit is seen to be: 


K = 5555.5 
If both tachometer and acceleration feedback are used the compensated 
system's characteristic equation becomes; 
st + 15S> + (150 + 5555.5 of ys? + (1000 + 5555.5 B YS + 5555.5 = 0 
where alpha = Ko beta = KR and K has been set at the stability limit. 


Parameter plane curves for this characteristic equation are shown in figure 


(3-23). 
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From these curves, the following analysis can be made. ‘The origin 
of the parameter plane corresponds to the roots of the uncompensated 
system. Since the zeta = O curve passes through the origin, two roots 
are located on the jw axis of the S-plane as was to be expected from the 
Routh array. The other two roots are also complex and are located at 
zeta = .8, and omega = 5. 

It is important to note that when an operating point involves two 
different pairs of complex roots, then the curves for two different values 
of omega and two different values of zeta will have to pass through the 
point. To determine which value of omega goes with which value of zeta, 
it is necessary to refer to the computer printout data. Due to lack of 
space, this data is not provided herewith. 

With K, set to zero, the effect of tachometer feedback alone corres- 
ponds to movement of the M point along the Y-axis. In figure (3-23) the 
unstable region is determined by an inspection of the way the constant 
zeta curves tend as zeta increases. Since the Y-axis is always in the 
unstable region, it is concluded that tachometer feedback alone cannot 
stabilize the system. 

The effect of acceleration feedback alone can be observed by moving 
along the X-axis of figure (3-23). If K, is varied between .01 and .06, 
the system will have two pairs of complex roots with the following ranges 
of values for zeta: | 

3 < zeta < .5 and .25 < zeta < .32 

If both tachometer and acceleration feedback are used, it is concluded 
that tachometer feedback will in general cause the zeta of one pair of roots 
to increase while the zeta of the other pair decreases. It appears that 


acceleration feedback alone would be better. 


From the graph it is determined that with K. = 0 and K .012, that 


complex roots are located at zeta = .45, omega = 4, and zeta = .32, omega 
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= 13. Since it is true that (.45) (4) = 1.8 << (2332) (3)82= 45a 
is apparent that the roots at zeta = .45 and omega = 4 are dominant. 


The specifications have therefore been met and the problem is solved. 


Example 3-19 
Problem: 

With the KjK, product set at the stability limit of the uncompensated 
system, use position feedback as shown in figure (3-24) to obtain domin- 
ant characteristic roots with maximum possible zeta and omega, and with a 
dominancy factor of about two to one (i.e., the ratio of the real part 
of any secondary root to the real part of the primary roots is about two 
or greater). 

Solution: 


From figure (3-24) the characteristic equation becomes: 


5° + 21008" + (X+ 1.2x10°)s + 100 of +@+10° = 0 


where alpha = KK, and beta = KK... Setting K, 


criteria, the value of KK, at the stability limit of the uncompensated 
system is found to be: 2.42%10°", 


= 0 and using the Routh 


In figure (3-25) are plotted parameter plane curves for the above sys- 
tem. Constant zeta-omega curves have also been included to assist in the 
dominancy considerations. 

The analysis proceeds as follows. The unstable region is to the 
left of the zeta = 0 curve. Root values for the uncompensated system 
are obtained for M point values along the Y-axis, since this corresponds 
to K, = Q. The stability limit of the uncompensated system corresponds 
to the intersection of the zeta = 0 curve and the Y-axis. The value of 
beta = KK, at this. point is observed by the Routh check. The M points of 
interest therefore lie along a horizontal line drawn through this point as 
shown in figure (3-25). 


The maximum zeta obtainable for any value of alpha along this line is 


about .33. As alpha increases then zeta-omega increases and the magnitude 
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of the real root decreases and the root becomes more dominant. Omega 
also increases with alpha. Therefore the overriding criteria is the 
dominancy factor. On this basis an M point located at alpha = KOK, 

= 2.5%10° and beta = KK, = 2.42 x 10” is chosen. Characteristic roots 
are read from the graph and are located at: zeta = .32, omega = 1680, 
zeta-omega = 550, and a real root at -1000. The dominancy factor is 
computed to be 550/1000 or 1.82. Since this is about two and the maxi- 
mum zeta and omega have been obtained, the solution is complete. 

It can be remarked that since KK, was set at a specific value for 
analysis purposes, the problem was then reduced to only one variable and 
root locus techniques could have been employed. This example illustrates 
the flexibility of the parameter plane and it also points out the interest- 
ing fact that root values determined along either a horizontal or a verti- 
cal line in the parameter plane constitute a root locus in terms of the 
variable parameter. Root values along a sloped straight line in the para- 


meter plane constitute a sort of nybrid root locus where the two root locus 


parameters are linear ly related. 


Example 3-20 
Problem: 

1. Set K such that the system of figure (3-26) will follow a 

ramp input, On = 1.0t. 

2. Steady state error should be less than two degrees. 

3. Step response overshoot should be less than 30%. 

4. Find values of § and P to meet the above specifications. 
Solution: 

It is known that steady state error = O,/K,» where in this problem 

=e K/10, and O, = 1.0. Error = 2° = .0349 rad. Therefore, Ke 28.6 


and K = 286. If the above system can be made second order dominant, the 


second order curves can be used. For an overshoot of 30%, a zeta of .35 
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is required. 


The characteristic equation of the compensated system then becomes: 


ag (11 ieee ne (10 + 1102 ce & (104 + 286 B )s + 286% =0 
In the above equation, alpha = p and beta = ¥ . Fourth order parameter 
plane curves are plotted in figure (3-27). It is seen that the zeta = 
.35 curve has a minimum value of beta = 255. Since beta is greater than 
ten, a multiple lead compensator or perhaps a lag-lead compensator is 
indicated. To enhance the chances of having to use no more than two 
sections, the minimum value of beta 25.5 is chosen for the single section. 
From figure (3-27) the following root locations can be read off correspond- 
ing to beta = 25.5 and alpha = 87: zeta = .35, omega = 8.2, and real roots 
at -90 and -4.5. Zeta-omega for the complex roots is three. Since 4.5 = 
1.5x3, the compensated system is second order dominant with a dominance 
factor of 1.5. 

At this point either a multiple section compensator can be designed 
from the above values of alpha and beta, or another scheme of compensation 


can be used. Both alternatives are considered in section (4). 
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4 Miscellaneous aspects of the parameter plane. 
4-1 Some general comments. 

In section (3-2-2) it was shown that constant zeta parameter plane 
curves of order two through five originate at a point where alpha = M/K 
and beta = N/K, where M, N, and K are determined by the zero and first 
power coefficients only. Inductive reasoning can be used to conclude that 
constant zeta curves of any order originate at this common point which is 
determined only by the zero and first power coefficients. An exception 
is when K = 0. In this case the origin point depends on higher order co- 
efficients and its location will be obvious given a specific equation. If 
K is not zero, the origin point is independent of the order of the charact- 
eristic equation. 

Inspection of the expressions for alpha and beta in section (3-2-2) 
indicates that the shape of the constant zeta curves as omega becomes 
larger is primarily determined by the coefficients of higher power, and 
in general the curves become more complex and less well behaved as the 
order of the characteristic equation increases. For a given character- 
istic equation, an increase in complexity can be observed as alpha and 
beta appear in more coefficients. As indicated in section (3-2-2) all 
constant zeta curves tend to plus or minus infinity. The relative magni- 
tudes of the coefficients determine whether the limit is plus or minus 
infinity. It is therefore necessary to choose a frequencyrange of inter- 
est before plotting the curves, thus limiting the analysis to one "window" 
of the infinite plane. 

Mitrovic curves, since they involve equations which have only one 
parameter appearing in each of two coefficients, are in general simpler 
and more well behaved than the parameter plane curves. It is also interest- 
ing to note that when parameter plane curves are plotted for characteris- 


tic equations of the Mitrovic type, the resulting curves are identical to 
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those plotted from the Mitrovic equations. This can be seen as follows. 
From reference (8) page 349, B and By are given by the following 


relations where notation has been changed to conform to this text. 


Y} 
k 
B aii Oa Bima. cee SON 
where B pa (2 4 r Bo) for ule? 
and 6, = 0, %, = ] 


Comparison of the tabulated @ functions in table (10-1) of reference (8) 


with appendix (IB) of this text leads to the following relation: 


k 
The By “B, Mitrovic characteristic equation is of the form: 
S" pee aGeCENBeS tape = 0 
2 1 oO 


If By = alpha and By = beta, alpha and beta can be computed from the 
parameter plane equations (2-10) and (2-11) and are as follows: 
alpha = Ss fea w¥ beta = S feupied k-ly 
K= 2 k k-l oe k k 


If the relationship between the 0. and U,. functions is employed in the 


Mitrovic expressions for BoB) one obtains: 


Boe 2a CD YY, B) = Zz dwt ¢-1) "0, 
It is seen that By = alpha and By = beta. The above procedure can be 
repeated for variables appearing in any two of the coefficients of the 
characteristic equation. This duality property is employed in section (4-2) 
to compensate parameter plane type characteristic equations employing normal- 
ized Mitrovic B By and B,-B. third order curves which have been plotted 
using the parameter plane computer program presented in section (6-1). 


Since in the parameter plane, a negative value for alpha or beta 


does not necessarily imply a negative coefficient of the characteristic 
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equation, one is not restricted to first quadrant values of alnha and 
beta, In most applications however a negative value for alpha or beta 
corresponds to a right half plane pole or zero in a cascade compensator 
or to positive feedback. For all of these cases only the first quadrant 
of the parameter plane is of interest. 

Since no stability criteria, either relative or absolute, has been 
established for the parameter plane, it is necessary to base the stability 
analysis on observing which way the curves tend as omega and zeta are 
varied. For this reason it is desirable to plot curves for as many values 
of zeta, omega, sigma, and if desired, zeta-omega, as is necessary to fix 


the pattern. 


4-2 Normalized third order curves. 
4-2-1. Discussion of the normalized curves. 

Third order curves are <vailable for finding the roots of polynomials 
where the coefficient of zero power is one variable and the coefficient of 
the first power is the other variable. Such a set of curves is presented 
in reference (8), and the curves are called normalized Mitrovic BBY curves. 

In this section a method is presented whereby parameter plane type 
characteristic equations of third order can be compensated on two different 
types of normalized curves. One type is the BBY curves of reference (8), 
and the other type is the BrB. curves which are presented in this section. 
The method here presented is divided into three cases. 

4-2-2 Derivation of the normalized transformations. 
Case I 


Characteristic equations of the type 


3 2 es / 
S + d,s Bye cle ¢ te Eee meet todo (4-1) 


are considered, where alpha and beta are the variable parameters. 


In equation (4-1) letting S = d,s one obtains: 
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ase —a 
d5s + dys + (b, te, Bt d))dos tbo of te fA d= 0 


The above equation is divided by fe resulting in: 


2 
3. 3 2 a 
s +s + s(b, 4 +c) Bit d,)/d, + (bo te. Q+ d /4, =Q (4-2) 


Let; 


2 


2 (4-3) 


Ba (b, o +c, e+ d,)/d 
3 
B (bof ter B+ d /d, 
Rearranging equations (4-3), two linear equations in alpha and beta are 
obtained: 


1 1 


bo + c, 6 =a qd 
e 


B. -d (4-3a) 
be+e B 
oO 


Ou 
NW KN NO 


- d 
O 9) 


Oo 
Equations (4-3a) can be solved by Cramer's rule or by the inversion of 


the matrix of coefficients to obtain: 


2 3 
alpha = ((d,B, - die. - (458, - de, /(,¢, - be) 


3 2, (4-4) 
beta = [(dB | - db, - (4, — db I/ (bc, - bc) 
Substituting equations (4-3) into equations (4-2) one obtains: 
ea = “i + B,Ss + B = 0 (4-5) 


Equation (4-5) is a normalized BB type equation. 


1 


In figure (4-1) parameter plane curves have been plotted for equation 


(4-5) where B, is the alpha or X-axis variable and Bis the beta or Y- 


1 
axis variable. 
Systems whose characteristic equation is of the case I type can there- 
fore be compensated on the normalized curves of figure (4-1), and the 
necessary values of alpha and beta can be found from equations (4-4). 
Case II. 
Characteristic equations are considered of the type: 


3 2 : 
S7 + (byt te, B+ d,)S" + (byt te, B+ 4))S+d,=0 (4-6) 
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» 3 * e @ e ® 
Letting S$ = do s in equation (4-6) and dividing through by d 
O 


one obtains: 


3 é , ai 2 3 es 3 ‘ ca 
s + (bh +c, G8+d,)a fd> + (boo +c, +d,)s/d? +1 = 0 (4-7) 


Let: 


3 


= ( 
B (b, A +c, iS + d,)/d. 


2 
, ¥3 Ca) 

8, = (b, +c, G + d)/d- 
Equations (4-8) when rearranged, reduce to two linear equations in alpha 


and beta as follows: 


it 
boo “+ Co oo dB. ° d., 


° (4-9) 
/ |= dpe e 
bio + Cy S dB, d, 
Solving for alpha and beta as in case I results in: 
5 
alpha = [(doB., - d,)c, - (dB, - de, I/ (boc, - b,c.) 
(4-10) 
3 $ 
= f a = = = 
beta i (dB, d,)b, (d 8. d,)b, 1/(b,c, bic, 
As a result of equations (4-8), equation (4-6) becomes: 
of + Bos. + B,Ss +1=0 (4-11) 


Equation (4-11) is a normalized B, 7B. type equation. In figure (4-2), 
parameter plane curves have been plotted for equation (4-11) where Be 
is the alnha or X-axis variable and By is the beta or Y-axis variable. 
Systems whose characteristic equation is of the case II type can 
therefore be compensated on the normalized curves of figure (4-2), and 
the necessary values of alpha and beta can be found from equations (4-10). 
Characteristic equations are considered of the type: 
53 4 (bof +e, @+4,)s" + (boo +e, Atd)Stbo Ate Ard =0 
(4-12) 


3 
® = q a e 2 2 
Letting S (b, + c, eG + d,)s and dividing through by (b, & + cy & “s d,) 
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equation (4-12) becomes: 


Me (biol +8 + 4,)8/(b,0% + CoP + ay)’ (4-13) 
ci 
+ (bo + cy B + d.)/(b, 0 + C4 + d,) = 


Let: 


wo 
ll 


(byl +e) B + 4,)/(b of + ey Bt ay) (4-14) 


B 
O 


3 
(bob +c) g + dd / (bh + cy 6 + d.) 


In equations (4-14) after rearranging, expanding, and collecting terms, 


one obtains: 


kz 
Byb5 of AE 1°78 + (2B aed, - b,)o4 : 2B b,c, A 8 if 
2 _ 
Bic d aaeCi) eemidoe” “1am 
and 
. 3 Be 3 2 
if 
Biba + B Aco 8 + id 22 *g + 6B oP 0° aa F + 
2 2 
3B .b6¢5 A + 3B bs do shi + 3Bo¢ an 6 + (4-15) 
(3B bd? - b Joh + (3B.c,45 - eB doe wel 20 
o 2 2 O o ZZ Oo o 2 O 


When the substitution of equation (4-14) is made in equation (4-12), a 
normalized By type equation results. Systems whose characteristic 
equation is of the case III type can be compensated on the curves given 
in figure (4-1), to obtain a value of By and By which can be substituted 
into equations (4-15). Equations (4-15) contains the solutions for alpha 
and beta. In the most general case where none of the coefficients in 
equation (4-15) are zero, a graphical solution of the second and third 
order polynomials can be made. That is the polynomials of equation (4-15) 
can be plotted on an alpha-beta plane and the necessary values of alpha 
and beta can be read from the curve intersections. When some of the co- 
efficients of equations (4-15) are zero, an analytical solution may be 


easiest, or perhaps a combination of the two can be used. 


4-2-3 Application of the method. 


Example 4-1 (Case l type) 
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Problem: 
In figure (4-3) find the values of K and K. to obtain: 
1. Characteristic roots at zeta = .5, and omega = 10. 
2. Error coefficient should be greater than or equal to 6. 
Solution: 
The appropriate characteristic equation is: 
5 4 1187 + (30 +o )S +6 = 0, where alpha = KK, beta = K. 
Specification 2. can be satisfied if: 
K. = K/G0+ KK) = 2 /CGQie) = > 
The characteristic equation is of the case I type, and the appropriate 
frequency transformation is therefore: S = d,s = lls 
Letting Wy represent the undamped natural frequency on the normalized plane, 
one can obtain via the frequency transformation: Wy = 91 


The value of zeta is unaffected by the transformation. From the 


curves of figure (4-1), a value of w. = .91 and zeta = .5 is seen 


2 N 


1 
to correspond to BY = .075 and By = .91. Employing equations (4-4) along 
with the appropriate coefficient values one can solve for alpha and beta. 
The result is: 


aiene =U(.91) 0 = 60 poceoechini ers eeiod 


pince alpha = KK, = 100K. = 80 it is seen that: K, = .8 and K = 100 = 


beta. 
The error coefficient is therefore: 
K, = f§ /(30 +L) = .91, which is greater ea .5 so the specifi- 
cations are satisfied. 
Note: If the error coefficient was required to be greater than .91, the 
specified value of zeta and omega would have to be modified so as to in- 


crease beta, decrease alpha, or both. 


Example 4-2 (Case II type) 
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Problem: 
In the system shown in figure (4-4) find K, Kee and K, to obtain 
the following: 
1. Characteristic roots at zeta = ./, omega = 10. 
2. The error coefficient should be greater than or equal to 6. 
Solution: 
The characteristic equation is: 
5? + (3 eo a + (2 + B )s + K = 0, where alpha = KK, and 
beta = KK, . The equation is of the case II type and the appropriate 
frequency transformation is: S = és, where d is as yet unknown. 
Therefore: 


I 


w= 10d (4-16) 


Employing equations (4-10) with the appropriate coefficients one obtains: 


\— 


WwW 


alpha =d By - 3 beta = 3, -2 (4-17) 


oO 


The error coefficient restriction can be satisfied if: 


3 
d= 12+ 68 or (d - 12)/6 = beta = cB, 2 (4-18) 


From equation (4-16) it is seen that; 
a = 10/w as = 100 /w2 and d_ = 1000/w> (4-19) 
Oo > “oO 7” Oo N 


Employing equations (4-19) and (4-18) one obtains: 


3 . 2 = 
1000/6w,, - 25 1000B, /w., - 2 or B, = 1.67/w, (4-20) 


1 


The B,-B, curves of figure (4-2) can now be employed. 
In figure (4-2) one can pick of various values of Wy along the line 
zeta = .7. The desired value of w, is one that satisfies equation (4-20). 


N 


Referring to figure (4-2), the following values can be obtained when zeta 


= ./: 
Wy = .8 and BS = 2.2 and using equation (4-20) it is found 
that By = 2.1. 
Also from figure (4-2) when zeta = .7 it can be seen that: 
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“ny = .76 and aa 2.2 and from equation (4-20) : Baus 222. 
The latter set of values are the desired ones and the corresponding value 
of Bo is 2.5. From equation (4-19): 
d_ = 1000/(.76)> = 2280 = K 
From equations (4-17): 
alpha = 29.9 beta = 379 
Therefore: beta = 379 = KK, = 2280K or K, = .166 and alpha = 29.9 = 
KK, = 2280K or K, = .0131. 
The error coefficient specification should automatically be satis- 
fied due to the way the analysis was performed, but as a check; 
K, = K/(2 + KK) = K/(2 +B ese 2 
Now: 
K/(2 + 8 ) = 2280/381 = 6 which confirms the result. 
Example 4-3 (Case III type) 
Problem: 
In figure (4-5) find values of gamma, P, and K to obtain: 
1. Characteristic roots at zeta = .7, omega = 15. 
2. The error coefficient should be equal to 20. 
Solution: 
The characteristic equation becomes: 
5? + (A+ 5)s* + (5% + 100 8 )S + 100% = 0, where beta = 
gamma, alpha = P, and K = 100 to satisfy specification 2. This is recog- 
nized as a case III type and the appropriate frequency transformation is: 
S = (< + 5)s. The appropriate curves are the BB, curves of figure (4-1). 
Employing equations (4-15) after dropping the terms with zero co- 


efficients one obtains: 


2 , : 
B, < “ + (10B, - 1) - 100 Q + 25B, = 0 (4-21a) 


B ee ise pao (75B_ - 100) + 125B_ = 0 
Oo O Oo Oo 


The latter equation when divided by B becomes: 
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om 341507 475 - 100/B )oL + 125 = 0 (4-21b) 


From the frequency transformation it is seen that: 
Wy = 15/(of + 5) (4-22) 
The procedure is to find values of Bo and Wa? employing the zeta = ./7 
curve of figure (4-1) such that equations (4-21b) and (4-22) are simultan- 
eously satisfied. This can best be done by trial and error. 
Trial l. 
From the curves, Bo = .065 and w= .58. From equation (4-22) one 
finds that alpha = 20.9. These values are used in equation (4-21b) to 
see if the left side of the equation becomes zero. Hence: 
(20.9)? + 15(20.9)* - 1465(20.9) + 125 = -14895# 0 
Teal 2. 
Try Ba = .07, then from the curves, Wy = -55. From equation (4-22) 
one finds that alpha = 26. Using equation (4-21b): 


(26)? + 15(26)* - 1355(26) + 125 = -17775 # 0 


Trial 3. 
Try BA = .075, then from the curves: Wy = .5. From equation (4-22) 


one finds: alpha = 25. From equation (4-21b): 
(25)? + 15(25)° =~ 1255(25) + 5(25) = -250, which can be consider- 
ed to be close enough to zero for the size numbersinvolved. Therefore 


from the curves, By = .45. 


Equation (4-2la) can now be employed to find beta: 


2 


beta = [(B, 4 ~ + (10B, - 1) + 25B 1/100 


1 


or 


[(.45)(25)* + (10x.45 - 1)25 + 25(.45)]/100 = 3.8 


beta 
The final results are: 


gamma = 3.8 P = 25 K = 100 
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4-3 Normalized parameter plane curves of higher order. 
In the preceding cection, essentially two types of transformations 
were made to normalize the third order characteristic equations. The 


first being the magnitude scaling, was tacitly assumed since the b. co- 


3 
efficient of the characteristic equation was taken to be one. Dividing 


through by b, constituted the magnitude scaling. The second transformation, 


3 
which was frequency scaling,was used to make one of the remaining coeffi- 
dents unity. This left only two variable parameters, and normalized 
curves could readily be plotted. If the characteristic equation is higher 
than third order, more than two parameters are involved since only two 
coefficients can be made unity, so general normalized curves of higher 
order than three are not feasible except in the special form as follows. 

In reference (5), Choe introduced normalized families of fourth order 
curves where the family parameter was taken as the B, coefficient. Para- 
meter plane transformations similar to those given in section (4-2) could 
be derived for the fourth order case, but they would be too complex to 


be of practical use. Obviously normalized curves for higher than fourth 


order are impractical, 
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4-4 Three dimensional parameter plane space. 
4-4-1 Discussion. 

Many compensation problems involve finding values for more than just 
two parameters. Such is true in multiple section cascade compensation and 
combination feedback and cascade compensation. This type of problem can 
be solved by conventional parameter plane methods if all but two para- 
meters are set at some arbitrary value. A more illuminating approach, 
however, involves the use of three dimensional parameter plane space. 
(When the problem involves or can be reduced to three parameters). 

Since the parameter plane equations are obtained by equating the 
real and imaginary parts of the characteristic equation to zero, unique 
solutions exist for only two parameters. However, a third parameter can 
be intorduced by plotting families of alpha-beta curves in two dimensional 
space. Theoretically a three dimensional parameter space surface could 
be plotted, but interpreting the results would be difficult. 

4-4-2 Example problem. 

In section (6-1) a computer program is presented which will plot 
families of parameter plane curves in terms of a third parameter as a 
variable. The third parameter may appear linearly or non-linearly in any 
of the coefficients. The following example shows one application of the 
method. 

Example 4-4, 

Experience has shown that if the parameter values necessary to com- 
pensate a system using a@ single section cascade compensator are not reali- 
zable, then the application of tachometer feedback will often permit the 
use of values within the acceptable limit. The system shown in figure 
(4-6) is the problem of example (3-20) but with tachometer feedback employed. 
Problem: 


Using the same criteria as in example (3-20), use tachometer feedback 
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to find more reasonable values for the cascade compensator parameters. 
Solution: 
From figure (4-6) it can be seen that: 


aan KP/(10P + KK, P) = K/(10 + KK, ) 


From example (3-20), K, = 28.6. 


Hence: K = 286 + 28.OKK, . 
The characteristic equation is found to be: 


sn (11 +% ys” a (10 + 110 + KK, B)S~ + (10¢ + KK, of + kK IS + KL = O 
where alpha = P and beta = gamma. The third parameter is seen to be KK: 
The above characteristic equation is seen to be identical to that of ex- 
ample (3-20) when KK. is set to zero. The error coefficient restriction 
can be incorporated into the above equation by letting K = 286 + 28.6KK, . 
This results in: 


5 + (11 +¢ ys? + (10 + 114% + + KK, BS” + (10o/+ KK, o + 2868 + 28. OKK P 


YS + 286% + 28 .6KK of = 0 
Third parameter values of 10, 50, 100, 200, 300, 400, 500, and 1000 
are investigated, employing the computer program of section (6-1). The 
resulting constant zeta and constant omega curves are shown in figure 
(4-7). For simplicity the only constant zeta curves that are plotted 
are for the required value of zeta of .35. If the M-point shown in 


figure (4-7) is chosen, the following values are applicable: 


KK, = 300 gamma = 10 P = 113 
This results in the compensator zero being located at S = -11.3, which is 
also a closed loop system zero. For the above parameter values the closed 
loop system poles were found using the digital computer and are as follows: 
-11.056 -90.835 -11.054 + 429.56 


The residue of the real root at S = -11.056 is approximately zero due to 


the closed loop zero at -11.3, so the complex roots which give the zeta 
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of .35 are dominant and the problem is solved. It is seen that tacho- 
meter feedback increases the undamped natural frequency of the complex 
roots by about a factor of three so the settling time is decreased by a 
factor of one-third. It is noted however, that settling time is not of 
interest in this example. 
At this point one should investigate to see if tachometer feedback 
alone could produce the desired results. The Routh check is first employed. 
The applicable characteristic equation is: 
s> + 118° + (10 + KK,)S + K = 0 
and 
K, = 28.6 = K/(10 + KK,); hence K = 286 ~ 28.6KK, . 
Using this value for K, the characteristic equation becomes: 
s? + 1187 + (10 + KK,)S + 286 + 28.6KK, = 0 


The Routh array is: 


1 10 KK. 

ll 286 + 28. OKK, 
-176 - 17. OKK, 0 
286 + 28. 6KK, 0 


Since a negative value of KK, is required to even stabilize the system , 
it is concluded that tachometer feedback alone will not work. 
4-5 Characteristic equations involving product terms of alpha and beta. 
4-5-1 Basic derivations. 

Assume that the coefficients of the characteristic equation are of the 
form: 


a = bote eth l@ +d, (4-23) 


Employing equations (2-7) one can obtain: 


0 (4-24a) 


eB, + @c, + “+P A, +D, 


of B, + Cee A, +D, = 0 (4-24b) 
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Where: 


mM 
> k, _k ae i aadle 
Ae GDh alee A, = Pa Gl) saaauy (4-25) 


UL. UL» Ba» C,> Dy» Bos Co» and D, are as defined in section (2). 
Equations (4-24) contain in general, two unique solutions for alpha 
and beta. Note that if A, - A, = 0, equations (4-24) reduce to equations 
(2-9) and determinants can be used to obtain the one unique solution. 
Solving equations (4-24a) and (4-24b) for alpha and equating the 
results one obtains: 
ete ie (4-26) 
Ca eee Tl eee a 


Equation (4-26) can be solved for beta resulting in: 


| 2 | 
betas 7 San + Age) + (Ayn t4aq)” - 4 Age A gy (4-27) 
——Ap_ Bo t= ADB AG BD 


: AC 
The deltas in equation (4-27) are shorthand notation for terms of the form: 
A gg = Byly - ByCys Sgn = A,Dy - ADD, 


Proceeding in a similar manner, alpha is found to be: 


| 2 | 
alpha = “CA, + SQ) = (Op, + Sg0) ' 4 Ag, 4, (4-28) 


24, 


, etc. 


Equations (4-27) and (4-28) contain four solution pairs for alpha and beta, 

only two of which satisfy equations (4-24). To find the two correct solue 
tions, equation (4-27) can be used to find two values of beta which can 
then be substituted into equation (4-26) to find the corresponding values 
of alpha. The two solutions can also be found from equations (4-27) and 
(4-28) where the two correct solution pairs are the ones that make equa- 

tions (4-24a) or (4-24b) go to zero. 
Constant zeta and constant omega curves can be plotted from equations 
(4-27) and (4-28). Proceeding as in section (2), conscant zeta-omega 


curves can also be plotted from equations (4-27) and (4-28) where By» Ch» 
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Dy» Bo Cos and dD, are defined by equations (2-15) and where: 
s a h 
al fe A, = Q (4-29) 
“ K=0 k*k-1 2 20 k “k 


Constant sigma curves however are no longer straight lines when an alpha- 


beta product is involved. In this case equation (2-17) becomes; 


; ee = k, gk Z k, ok 
AE (D6 5+ PZ CDi 66+ KP Z (1h S © 4-30) 
Se k k 
+ & (-1)"d. gs = 0. 
K=0 


If one assumes the following notation: 


mM Ww 
pp = 2 Ena c* 6cClae tebe 
K=O k K=0 k 
mM yi 
BBB= = ab ae BO= = (aan ee 
K=0 k k=0 k 
Then equation (4-30) becomes: 
of BBB + 6 CCC + & & BC + DDD = 0 (4-31) 


Equation (4-31) is a special form of a conic section and can be plotted 
by solving for alpha and incrementing beta over a range of values of interest. 
The above equations are programmed for the digital computer in section (6-2). 

Example problems involving the above concepts are presented in the 


following section. 





4-6 Design of double section cascade compensators. 
4-6-1 Discussion. 

Double section compensators can be designed by use of the parameter 
plane in two ways. The double section compensator can be made equivilant 
to a single section at a specific complex frequency of interest. Here 
the system is first designed using a single section compensator but the 
single section parameter values turn out to be physically unrealizable. 
This method has the disadvantage that control is maintained over only one 
pair of complex roots and dominance is difficult to ensure. 

The second method involves writing the characteristic equation with 
a double section compensator inserted and drawing the parameter plane 
curves. Control over all the roots is obtainable and the dominance pro- 
blem is much simplified. 

Both the above methods involve characteristic equations with alpha- 
beta product terms appearing in the coefficients, and new parameter plane 
equations have been derived to handle this situation in section (4-5-1). 
4-6-2 Design of a double section compensator on the basis of given single 

section parameter values. 

This method was proposed by Hyon in reference (6) but Hyon used the 
Mitrovic equations to obtain a solution. Parameter plane techniques will 
now be employed. 

Let the open loop transfer function of a double section cascade 
compensated system be: 


Ya Sug, St Pal i) 5 ES eel 5 (4-32) 
(S + PCS + Po) 


The uncompensated system's forward path transfer function is G, 
where unity feedback is assumed. The open loop transfer function of a 


single section compensated system is: 


(S++ P/ Seow 1Gaeeel (4-33) 


(S + P) 
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Assume conplex roots are required at: 


| 2 
Sy = - 7 4) + jw, lL - Va (4-34) 


Equations (4-32) and (4-33) are equated and G is divided out: 


Ee (S+Pi/ ¥ )GS+ PSX 4) (S + P/X ) (4-35) 
“(S + P (Ss - ey (S + P) 


After rearranging, collecting terms of like power, and dividing by S, 


equation (4-35) becomes: 


(4-36) 


P, ¥ 4] + Pot SHRP, tePjoceBaeS ooneBQpe- P,P, = 0 
The unknowns in equation (4-36) are Py Po» ¥ 1? Yo and S. 
By substituting a specific frequency for S and equating the real and imagin- 
ary parts of equation (4-36) to zero separately, two equations in two un- 
knowns can be obtained. The parameters 4 and X 5 can be pre-set to a 
fixed value, thus determining whether the compensator will be double lead, 
double lag, or lag lead. Values for the remaining unknows Py and P, can 
then be obtained from the resulting two equations. 

It can be noted however, that when 84 and ts are given fixed values, 
the parameter plane techniques of section (4-5) can be applied directly 
if one lets P, =f , ES p , and P,P, =o@ . Equations (4-27) and 


(4-28) give the desired values of Py and P, when the specified value of 
zeta and oneas for the complex roots is substituted. Only one or at most 
two points in the parameter plane are of interest thus obviating the need 
for plotting curves. Hence straight analytical techniques can be used. 
Example 4-5. 


Problem: 


Apply the above techniques to design a double section filter equivalent 
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to the single section filter of example (3-20) at the specified frequency 
of zeta = .35 and omega = 8.2. Of the possible solutions, choose the one 
that makes the specified complex roots most dominant. 
solution: 
From example (3-20) a value of gamma = 25.5 was needed, indicating a 
double section lead filter might work. The value of P was 87. Let % a 
85 = 5, which is a reasonable value for a lead filter. When the above 


values are substituted, equation (4.36) becomes: 


s* 4 (210 + 21 B- 2118)8 + 2% - 348( 4+ B) = 0 (4-37) 


when zeta = .35: U,=-1; U. = 0; U, = 13; U, = .7 
-1 oO 1 


From equation (4-37) the coefficients can be used along with equations 


(2-10) and (4-25) to obtain the following quantities: 


A, = -25 A, = 0 
BS = 348 By = -172 
Cc, = 348 Cc, = -172 
D, = 67.2 Do = 17397 


The deltas are then found to be: 


D gp = 4.349 x 10° A ne = ~6.066 x 10° 
A ac ="0 DA ny = 4-349 x 10° 
A ae = 4300 A, = ~4300 

A. = 6,066 x 10° 


Using equations (4-27) and (4-28) the solutions are: 


x, = 16.7 of 4 = 84.4 6, = 84.4 B . = 16.7 
From equations (4-24) it is found that 04 7 paired with fs 1 and 0% 9 
paired with & 9 are the consistant solutions. 
Since in equation (4-37), the coefficients of alpha and beta are 
identical, both solutions produce the same result. 1 and 6 1 are arbi- 


trarily chosen. The characteristic equation of the compensated system 
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then becomes: 


3 42 5 


) TY 5.9678) % 10rS5 Se Teen ae 


S” + 112.148" + 2533.28 
x 10° = 0 


Using the digital computer the roots are found to be: 


-4.205 -16.79 -85.51 -2.82 + j7.674 
The complex roots are at the specified value of zeta and omega. 

Comparing the above roots with those obtained in example (4-4) where 
the single section plus tachometer feedback was employed, one concludes that 
the double section compensator produces dominant complex roots whereas the 
tachometer feedback and single section scheme does not. In example (4-4) 
however, the complex roots were found to be effectively dominant since the 
residue of the nearby real root was about zero. In this example, closed 
loop zeros are located at -3.34 and -16.9. Therefore, the magnitude of 
the residues of the complex poles and the real pole at -4.205 are almost 


the same so it is fortunate that the complex roots are dominant. 


4-6-3 Design of a double section compensator using general parameter 
plane methods. 

This method involves incorporating the double section cascade compen- 
sator equations into the uncompensated system's characteristic equation. 
This technique has the advantage that parameter plane curves can then be 
drawn and control i: maintained over all the characteristic roots rather 
than only the two specified complex roots as is the case with the method 
of section (4-6-2). 

Here the technique can best be explained by an example. 

Example 4-6 
Problem: 

Figure (4-8) shows the system of example (3-20) but with a double 

section compensator employed. Solve example (3-20) using a double section 


compensator as indicated in figure (4-8). 
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Solution: 
For comparison with the results of example (4-5) let: 
bo ae ee 
The characteristic equation then becomes: 
so 4 (ht A+ yst+ ( AB + 11o +11 + 10)” + 
(idG +10% +1084 7150)S° + (LOX B + 14300 + 
1430 @ )S+ 2608 =0 (4-38) 


Here; OO =P r moe and AB = P 


ql” Io 

Parameter plane curves for equation (4-38) have been plotted in 
figures (4-9) and (4-10) utilizing the computer program of section (6-2). 
Since the curves are rather complex, the constant zeta and constant omega 
curves are plotted in figure (4-9) and the constant sigma curves are plot- 
ted in figure (4-10). In figure (4-9) it is noted that there are two 
sections of constant zeta curves for each value of zeta and two section of 
constant omega curves for each value of omega. This agrees with the re- 
sults proven earlier that two unique solutions exist for each value of 
zeta and omega. 

In figure (4-10), the discontinuities of the constant sigma curves 
are indicated by the horizontal straight line. The straight line is not 
part of the curve and should be ignored. A study of the curves indicates 
that for a given curve, if one chooses a point, say (o4 © GE D>»? then 
there exists a corresponding point ( B ? of Dp: This is due to the 
fact that the coefficients of alpha and beta are identical in the 
characteristic equation. 

If the M-point indicated in figures (4-9) and (4-10) is chosen then 
the following characteristic roots are indicated: 

-9.21 | -~26.14 -65.34 -2e1 + 35061 

The accuracy of the root values was obtained from the printed out 


computer data, given that alpha = 31.2 and beta = 62.5. The complex roots 


are located at zeta = .35 and omega = 6. 
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The roots obtained in example (4-5) are: 

-4.205 -16.79 -85.51 -2.82 + 47.674 

Comparison of the two sets of roots indicates that in example 
(4-5) the ratio of the nearest real root to the real part of the complex 
roots is 1.44 whereas in this example it is 4.4. Hence the dominancy 
factor has improved by a factor of three in this example. 


The effectiveness of the method of this section is thus apparent. 
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5 Root locus digital computer programs. 

The programs of this section were written using Fortran 60 along with 
subroutines available at the computer facility of the U. S. Naval Post- 
graduate School. It is necessary to have the coefficients of the poly- 
nomial or characteristic equation to use the root locus program of section 
(5-2), so a program is presented in section (5-1) to compute the coeffi- 
cients in case the equation is in factored form. 

The use of the programs is explained in the comment cards at the 


beginning of the programs. 
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6 Parameter plane digital computer programs. 

The programs of this section were written using Fortran 60 along 
with subroutines available at the computer facility of the U. S. Naval 
Postgraduate School. 

The use of the programs is explained in the comment cards at the 


beginning of the programs. 
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/- The complementary roles of the parameter plane and root locus. 

As was mentioned previously, parameter plane curves are infinite in 
extent. The curves also exhibit discontinuities and become less well 
behaved as the order increases. In addition, choosing a good graph scale 
involves a good deal of trial and error or pre-plotting calculations. 

The root locus can be very useful in overcoming some of the above diffi- 
culties. 

The general procedure is to use the root locus to limit the range 
of parameter values of interest, and to gain additional insight into the 
problem. This is particularly true if a computer is used to plot the root 
locus. Since a correct root locus graph scale can be fixed in one or two 
computer runs, several root loci can be plotted in fairly short time. 

The parameter plane can then be used to complete the problem. 

The technique can best be illustrated by the following example. 
Example 7-1 
Problem: 

For the system shown in figure (7-1) find values for Ky: Ky» Ky» KR 
and K. to give a good transient response. Low settling time is the 


primary consideration. Suggest modifications to improve the system. 


The following parameters have fixed values: 


P, = 6.28 x 10° Py = 2.5 x 10° oa °68 eee de 
zeta = 0.5 
Solution: 


‘It is known that the block diagram elements of figure (7-1) that 
involve Ky and Ke were added on the basis of physical reasoning to help 
reduce the settling time. The basic uncompensated system consists of the 


single loop unity feedback control system involving Ky» Ky» and K,. From 


the block diagram reduction of figure (7-1) the characteristic equation 
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is computed and the given parameter values are substituted. Letting 


alpha = K hy and beta = KKK, one then obtains: 


5 7.4 13 20 


5 + 1.281 x 10S + (3,194 xe . oye +7 (2,0/55xenu 


25 


+O + 2.5 x 107 of. sg” + (6.18 x 10°” + 6.28 x 10°B a 


, | 13, _ 
K,KK,K,K.)s + 3.94 x 10 6 = 0 (7-1) 


To study the effect of the Ky compensator, Ke 


in equation (7-1) and parameter plane curves are plotted in figure (7-2) 


is set equal to zero 


for the remaining system. 

The unstable region for these curves is above and to the right of 
the zeta equals zero curve. From the Routh check, the stability limit 
gain of the uncompensated system, i.e., K KOK,» is equal to 1.025 x 107°. 
This value corresponds to point B in figure (7-2). The effect of the Ky 
compensator is to make the system less stable as KK, is increased. With 
KKK, set to 1.025 x fo. a root locus is plotted in figure (7-3) with 


KK, as the variable. This plot shows the system to be unstable for all 
KK, greater than zero. Gain changes have small effect on the right half 
plane root locations. 

Now KK, is arbitrarily set to 9.1 x 10° and KKK, is left unchanged. 


The effect of the K, compensator is seen from figure (7-4) to have a 
stabilizing effect. It is known that settling time is inversely propor- 
tional to the undamped natural frequency of the complex roots. In figure 
(7-4), point A represents a zeta and omega that would be satisfactory for 
the complex roots, but the location of the real root shows the latter to 
be dominant. 

Two avenues of approach appear to be open. Try a different means 
of compensation, or modify the suggested means. The former is investi- 
gated first. 


The effect of tachometer feedback around the entire forward path, 
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with Ky = K, = 0, is shown in figure (7-5). From the root locations as 
tachometer gain increases from point A to B to C, it is seen that the 
real root is highly dominant. Tachometer feedback would therefore be un- 
satisfactory. Since the system's operating frequency is in the mega- 
cycle range, higher forms of derivative feedback would not be practical. 

The fact that the Ky compensator, which is the feedback path around 
KK, in figure (7-1), makes the system less stable, suggests that this 
path could be eliminated. When this is done the resulting system is 
shown in figure (7-6) and the characteristic equation is as follows: 

5 7.4 13 


Seto. 2a x 10 S$ + 82 1945x810 s? + (26.67 x 10°7+ 8 Wea te 


25 


(6.18 x 10°? + 6.28 x 10°B +oL )s 3.94 x 1098 =0 (722) 


In equation (7-2), alpha = Kj KoK3K, Ke and beta = K, KK, - 

The advantage of this scheme over the original system is that the 
KK compensator can be realized by an R-L-C circuit, whereas due to the 
physical nature of the problem, the original scheme would require electro- 
mechanical implementation with the inherent disadvantages. 

In figure (7-6), KKK is again set at the stability limit of the 


uncompensated system and a root locus with variable a is plotted in 
figure (7-7). A study of figure (7-7) shows that a slightly better domin- 
ance factor can be obtained with the modified system. It is concluded 
that the latter system not only performs as good or better than the former, 
but it is simpler and cheaper to implement. 

To see if system performance can be improved by increased gain, the 


forward path gain K is increased by a factor of ten. The resulting 


Po 
root locus in figure (7-8) indicates that the system is unstable for all 
K, Ke This illustrates the difficulty in choosing the best values for 

system parameters by root locus techniques when more than one parameter 


is involved. 


The parameter plane can now be advantageously employed. The resulting 
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curves are plotted in figure (7-9) with variables K KOK, and K,KoK3K)K. 
A close study of the curves indicates that the M-point shown is perhaps 


the best one. The five characteristic roots can be read directly from 


the curves and are as follows: 


zeta oO omega = 3 x 10 


zeta 


1} 


edu omega = 6.5 x 10 


S = -2.2 x 10° 

The damping factor of the first pair of roots is 1.65 x 10° which 
when compared to 2.2 x 10°, indicates that these roots are dominant. 

In the original uncompensated third order system, for a zeta of .55, 
the maximum obtainable omega for the dominant roots is only about .3 x 
10. On the basis of the dominant roots only, one can conclude that the 
compensation has reduced the settling time by a factor of ten. The trans- 
ient response of the compensated system to a unit step input was computed 
by digital computer and is given in figure (7-10). The settling time is 


about four microseconds and the overshoot is 50%. 
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8 - Conclusions. 

Parameter plane techniques have been applied to the compensation of 
linear feedback control systems. In particular, general equations have 
been derived for the cases of feedback, cascade, and combination feedback- 
cascade compensation, to enable one to analytically place a pair of com- 
plex roots at a specific location in the S-plane, while simultaneously 
satisfying the steady state accuracy requirements. 

A dominancy technique has been introduced whereby once the specified 
pair of complex roots is fixed, the remaining roots of the characteristic 
equation can be placed to ensure that the specified roots are dominant. 

Sketching techniques are developed enabling one to quickly sketch the 
zeta equals zero and zeta equals one-half curves. This can be useful in 
determining the type of compensation to use and in choosing an appropriate 
graph scale for the digital computer plots. 

Graphical solutions on the parameter plane are discussed in terms 
of engineering example problems. 

Miscellaneous aspects of the parameter plane are discussed. In parti- 
cular, transformations are derived to enable one to compensate parameter 
plane type characteristic equations on normalized Mitrovic third order 


B38 and B curves. Three dimensional parameter plane space is dis- 


1 1779 
cussed and applied to an example problem. 

A derivation of parameter plane equations involving product terms of 
alpha and beta is made and the results are applied to the design of double 
section cascade compensators. Double section compensators are designed on 
the basis of unrealizable single section parameters and by incorporating the 
double section compensator into the characteristic equation and plotting the 
parameter plane curves. 


Digital computer programs are introduced that one can use to plot root 


loci and parameter plane curves. Parameter plane curves can be plotted for 
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characteristic equations involving three parameters and alpha-beta pro- 
duct terms. 

Finally an engineering example is presented that points out the 
complementary nature of the root locus and parameter plane. 

A basis for further investigation involves plotting constant band- 
width curves on the parameter plane and determining the nature of para- 
meter plane curves resulting from characteristic equations containing 


squared terms of alpha and beta. 
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APPENDIX II 


A. Synthesis of R-C lead network 


Ci 
Lin R R= Ro 
Fig. a 1 
In figure (a), = = — St ie = 7 
IN S+ez,[ - | See 
where W2R,¢C, and olf Re /( R,+ Re |) . The above transfer function 


= 


has aD. C. gain of of so to make the D. C. gain unity an amplifier of Gy 
gain will have to be added. The pole to zero ratio, Ly is greater than 
one, indicating that the circuit of figure (a) is a lead network. 
If the pole to zero ratio and “[ become very large the filter be- 
haves like a pure differentiator and noise problems arise. If the pole to zero y 


ratio s kept less than ten, the noise problem is reduced. 


B. Synthesis of R-C lag network 


R, 9 
Lin ; Lo 
ee eee 
Fre. 0b 


In figure (b) the transfer function 1S: 


R | 
S| Bs os Ce( hit a i eae + ae 
ea Ss 
s C2(RitRz) He 


where T= Ca(Ri+ Pe) and of = Re / (Rite). 


The D. C. gain of the above transfer function is unity and the pole to 
zero ratio is alpha. Since alpha is less than one, the circuit of figure 


(b) is a lag network. 


75 





BIBLIOGRAPHY 


Ross, E. R., Warren, T. C., and Thaler, G. J., Compensation Based 
on the Root Locus Approach, AIEE Paper, Sept. 1960. 


Pollak, C. D., and Thaler, G. J., S-Plane Design of Compensators 
for Feedback Systems, IRE Transactions, May 1961. 


Mitrovic, D., Graphical Analysis and Synthesis of Feedback Control 
Systems, IEEE, Theory and Analysis, II Synthesis. AIEE Trans- 
actions, Part II, Jan. 1959. 


Ohta, T., Mitrovics Methed - Some Fundamental Techniques, Research 
Paper No. 39, U. S. Naval Postgraduate School, Jan. 1964, 


Choe, H. H., Some Extensions of Mitrovic's Method in Analysis and 
Design of Feedback Control Systems, M. S. Thesis, U. S. Naval Post- 
graduate School, 1964. 


Hyon, C. H., A Direct Method of Designing Linear Feedback Control 
Systems by Applying Mitrovic's Method, M. S. Thesis, U. S. Naval 
Postgraduate School, 1964. 


Siljak, D. D., Analysis and Synthesis of Feedback Control Systems in 
the Parameter Plane, Part I, Linear Continuous Systems, IEEE Trans- 
actions, Nov. 4, 1964, 


Thaler, G. J., and Brown, R. G., Analysis and Design of Feedback 
Control Systems, McGraw-Hill Book Co., Inc., 1960. 


176 


























— 7 — 



























eh ; 
a eB a 
ae Pv AE 
> yh ai (At Mere 
LY ory a Felt iy on NE Pet yt ‘a f 
; ; Ky : Aa y (ht 3 1. EP uw r! 
ek it LP ire ee 5 
ve Mit i Ny WD ea f na 1 a ee : 
= 5 | vi Ai bah , =} Pa . \ L, 
4 ry he ioe 4 i arr 
ae . fl y : ey 4) N i i] 7, 
AA a ‘ ; a7 ne = mp” 
H tt Fy vi i, u F 
1 ) ' i 1 1 " 
arin ae vee 4 on 
* i Uy Th | , i iy 
rN te gt Fae Ne wire ' 
t f oa i ‘| dy P Peay 
i My ae ve 7 i, 
i] ¢ i 
AM ‘FD ie ‘ . 
i i r J 
t 77 ' i af A 
> i ; if { F , 1 
a a i ' 
i a ve ! 
~ t= . i oar Pe f 
~ ro —» - i] | 1 4. 7 i : 
rs me au , _ 4 - ‘ ou Ma ' | Af 
‘ 1% ee - J “Se A , ' aM, 
t- 7 — Fe AA abe ' 
t | a i} 
a a t wT) yy 
= - Ao “ “wa ny ‘ y i 
~ =~ as A - i rl i , f 
=_ = - “= oy ' 17 , ‘ 
_ ; 1 i . i R 
_ - i} | ' | a 
-— - ea 
= \ iD 
T ass | | - ' 
| Ae a | i} 
is 
= - 1 
| 
ir \ 
| 
. i 
4 j 
VY 
i] i] 
V 
‘ 
Ta 
/ " 
L # 
i 
et 
te 5 ' 
J L 
i 
i; U i 
i ie PF ; ' , 
' 
y Ti i) "i 
| ' 
' | | 
i i F ( 
i 
te Le i 
ad 1 \ ) 
Poy 1% ‘ , 
Nes 7% 
TL o U i \ if ! ' 
t i} 1 =i) Fh ‘ et A yal : 
M ] : 
¥ a Vel ' H 0. . 
r , i i 
_s Aye "y H 
7 eh ve i} i} 
; i , aM a7 Pw | ' a + 
ao 1 ; re | 1 BF | A r) a? te - v 
fhe oa af ; H Ta 1 i 
A y Wir | 7 at Pere 
e Ly i AS # if - vt | } 
‘ ' a 
i ak 4 Ve f ‘ 
I tee f . 
ee te UE ~~ io 10 
er a 1, 
“ i a ve Th pee nn ye 
7) att = oiak ff faba Me oe Pa as “a COE) AD A Ye a 
iv Lf 1 eae VAL wi ’ i Pie ani Ue ey Wh 
i i . hy EE ol . - i veh , Th 
a re yeas Po tek: n doae 7 AA! 
Mf ’ Ta yp a 


__ 
ae 7 ~at, 7 







Ay ee 


eae 


eee 
tie 
a 





i 


ive an ee 


7 


a f L | yl E 
tek. Wis ee) t 1} , 
ce fy 7 a eX rh iy a by aie a . rr eS 4! 4 


i. ie i ee | 


